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ABSTRACT 

We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf 
spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show 
significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in 
Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the 
Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon 
stars. However, the more metal-poor stars in the current sample appear to be generating less dust. 
These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC 
dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and 
infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar 
to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study 
of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient 
quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star's 
lifetime on the asymptotic giant branch. 

Subject headings: circumstellar matter — infrared: stars — Magellanic Clouds — Local Group 



1. INTRODUCTION 

Stars on the asymptotic giant branch (AGB) are an 
important source of dust injected into the interstellar 
medium in the Milky Way (e.g., Gehrz 1989; Habing 
1996). How important they are in more metal-poor en- 
vironments is an open question, with consequences for 
the early history of the Milky Way, current conditions in 
other smaller Local Group galaxies, and even for galaxies 
in the high-redshift Universe. 

The sensitivity of the Infrared Spectrograph (IRS; 
Houck et al. 2004) on the Spitzer Space Telescope 
(Werner et al. 2004) has made it possible to explore this 
question by observing the dust forming around individ- 
ual evolved stars in environments spanning a range of 
metallicities, both in our own galaxy and elsewhere in 
the Local Group. The IRS has observed dust around 
AGB stars in the Large Magellanic Cloud (LMC; Zijlstra 
et al. 2006; Buchanan et al. 2006; Leisenring et al. 2008), 
Small Magellanic Cloud (SMC; Sloan et al. 2006; La- 
gadec et al. 2007), Fornax Dwarf Spheroidal (Matsuura 
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et al. 2007), Sagittarius Dwarf Spheroidal (Lagadec et al. 
2009)10, Sculptor Dwarf Spheroidal (Sloan et al. 2009), 
and several Galactic globular clusters (Lebzelter et al. 
2006; Sloan et al. 2010; McDonald et al. 2011). These 
metalliticies of these systems span a range of —2.1 < 
[Fe/H] < 0.0. Comparisons of these samples with each 
other and with samples from the Galactic disk reveal 
that the amount of dust produced around oxygen-rich 
AGB stars decreases in more metal-poor environments 
(Sloan et al. 2008, 2010), but for carbon-rich AGB stars, 
the amount of dust observed shows no measurable depen- 
dence on metallicity (Matsuura et al. 2007; Groenewegen 
et al. 2007; Sloan et al. 2008). 

Stars on the AGB will produce either carbon- 
dominated or oxygen-dominated dust, depending on 
their dredge-up histories and initial abundances. AGB 
stars generate carbon via the triple-alpha sequence 
(Salpeter 1952) in a helium-burning shell around an inert 
C/0 core (e.g., Iben & Renzini 1983). The helium fusion 
proceeds in a series of thermal pulses in the interior, lead- 
ing to pulses of convection which dredge newly produced 
carbon up to the surface. The dredge-ups can raise the 
photospheric C/0 ratio above unity, converting an AGB 
star with a spectral class of M giant to a carbon star. The 
envelopes of AGB stars are unstable to pulsations with 
typical periods of hundreds of days, making them read- 
ily identifiable as long-period variables (LPVs). These 
pulsations in the envelopes may even drive the mass-loss 
process (see Mattsson et al. 2008, and references therein). 

The formation of CO in the resulting outflows will ex- 
haust all of the available carbon or oxygen, whichever is 
less abundant, leading to a chemical dichotomy in the 
dust which will condense out of the outflowing gas. Alu- 
mina and silicates will dominate the shells around M 

Also known as the Sagittarius Dwarf Elliptical Galaxy, or 
SAGDEG, and not to be confused with the Sagittarius Dwarf Ir- 
regular Galaxy, or SAGDIG. 
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giants, and amorphous carbon will dominate the shells 
around carbon stars (e.g., Martin & Rogers 1987; Onaka 
et al. 1989; Egan & Sloan 2001). 

In more metal-poor galaxies, stars of lower initial mass 
will become carbon stars on the AGB. Counts of carbon 
stars in the LMC and SMC reveal this fact observation- 
ally (Blanco ct al. 1978, 1980; Cioni & Habing 2003), and 
it is expected theoretically (Renzini & Voli 1981; Karakas 
& Lattanzio 2007). A recent infrared census of the SMC 
with Spitzer reveals the consequence: carbon stars pro- 
duce more dust than their oxygen-rich AGB counterparts 
or red supergiants (Matsuura 2012), possibly much more 
(Boyer et al. 2012).^^ The SMC serves as a proxy for 
metal-poor galaxies too distant for their constituent stars 
to be studied individually. By studying even more metal- 
poor galaxies in the Local Group, we can push to even 
more primitive systems. 

The infrared spectra of seven carbon stars beyond the 
Magellanic Clouds have been published so far, six in For- 
nax (Matsuura et al. 2007), and one in Sculptor (Sloan 
et al. 2009). This paper presents spectra from the IRS 
for a larger sample of carbon stars in these two galax- 
ies, as well as the Carina and Leo I dwarf spheroidal 
(dSph) galaxies. The stars in the sample arc quite faint 
in the infrared, and the development of a new algorithm 
to extract spectra from the two-dimensional IRS images 
(Lebouteiller et al. 2010) makes it possible to analyze 
their spectra similarly to closer and brighter samples. 
Additionally, near-infrared (NIR) monitoring from the 
South African Astronomical Observatory (SAAO) has 
provided information on the pulsation modes and peri- 
ods of the targeted stars which was not available when 
Spitzer observed them. These two improvements give us 
the opportunity to extend the previous comparisons of 
mass loss and dust production in evolved stars to more 
distant galaxies with lower mctallicities. 

Section 2 presents our targeted galaxies, with an em- 
phasis on their distances and metallicities, and explains 
how we selected our sample of stars. Section 3 describes 
the observations and data reduction. In Section 4, we de- 
termine bolometric magnitudes and use these to re-assess 
the metallicities. Section 5 presents the spectroscopic re- 
sults, and Section 6 discusses the consequences of our 
findings. 

2. THE SAMPLE 

2.1. Target galaxies 

Our targets sample the evolved stellar population in 
four dwarf spheroidal galaxies in the Local Group. Sculp- 
tor was the first dwarf spheroidal discovered, quickly fol- 
lowed by Fornax (Shaplcy 1938). Leo I was uncovered 
during the first Palomar Sky Survey (Harrington & Wil- 
son 1950), and Cannon et al. (1977) detected the Ca- 
rina dwarf while conducting the Southern Sky Survey 
from the European Southern Observatory (ESO). Table 1 
presents some basic data for these galaxies that we will 
use for the remainder of the paper, and it requires some 
explanation. 

The distance moduli were determined with weighted 

The contribution from SNe at this time is highly uncertain 
due to contradictory measurements at different wavelengths; see 
(Matsuura et al. 2011) for the possibility that SNe can produce 
large amounts of dust. 



averages of published distances based on standard can- 
dles such as RR Lyrac variables, the horizontal branch 
(HB), and the tip of the red giant branch (RGB). The 
uncertainties in the distance moduli are statistical and 
do not reflect systematic errors, which are likely to be 
larger. Table 2 lists the individual distance measure- 
ments used to determine the results in Table 1. The 
entries in Table 2 are not meant to be exhaustive; some 
measurements rendered redundant or obsolete by more 
recent work are not included. Uncertainties smaller than 
0.04 magnitudes have been raised to that limit for the 
purpose of weighting (and noted with parentheses in Ta- 
ble 2). 

Table 1 includes our assumed values of interstellar 
reddening from foreground extinction in the Galaxy, 
E{B — V). In this study, they only influence our de- 
rived bolometric magnitudes, and the influence is small, 
because the extinction is small to begin with and smaller 
still in the NIR, where the carbon stars emit most of their 
energy. As an example, a reddening of 0.03 magnitudes 
corresponds to an extinction of 0.03 magnitudes at J and 
0.01 at K (using the extinction law of Rieke & Lebofsky 
1985), making the impact of reddening smaller than our 
uncertainty in distance. Our assumed reddening values 
are consistent with the literature and the infrared dust 
maps of Schlegcl et al. (1998), except that reddening de- 
rived from the dust maps is higher than typical values 
used for Carina (0.06 vs. 0.025 mag.; e.g., Smecker-Hane 
et al. 1994; Mateo et al. 1998). 

2.2. Metallicities 

Our objective is to understand how the infrared spec- 
tral characteristics of carbon stars vary with metallicity, 

making it important that we understand the metallicity 
distribution functions (MDFs) of the parent populations 
of our targeted carbon stars in each galaxy. Tackling 
this question is not easy, given the complex star forma- 
tion histories of these systems. Each galaxy is unique in 
this regard. 

2.2.1. Sculptor 

Sculptor has two populations with different ages, 
metallicities, and spatial distributions. Its color- 
magnitude diagram (CMD) shows two RGB bumps and 
two HBs, consistent with two populations with distinct 
metallicities (Majewski et al. 1999). Hurley-Keller ct al. 
(1999) noticed that the red HB, which arises from the 
more metal-rich population, is confined to the center of 
the galaxy. 

In the past decade, multiple studies have used the Fi- 
bre Large Multi-Element Spectrograph (FLAMES) at 
the Very Large Telescope (VLT), measuring the Ca II 
triplet (CaT) at 0.85 /im to determine the metallicity of 
hundreds of RGB stars. Those that focus on the core 
regions (within ^10-15' of the center^^) tend to pro- 
duce more metal-rich MDFs compared to those sampling 
larger parts of the galaxy. Inside 12' of the center, Tol- 
stoy ct al. (2004) find a MDF with <[Fe/H]> = -1.49 
± 0.35, but outside that ellipsoidal radius, <[Fe/H]> = 
—1.91 ± 0.27 (our calculations based on figures in their 

All distances in this section are ellipsoidal, or foreshortened 
away from the major axis. See Irwin & Hatzidimitriou (1995) for 
geometrical details. 
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TABLE 1 

Dwarf Spheroidal Galaxies Studied 



dSph 


Distance 


Adopted 


Mean 


Adopted 


galaxy 


modulus 


E{B~V) 


[Fe/H] 


[Fe/H] 


Sculptor 


19.64 ± 0.04 


0.02 


-1.56 ± 0.40 


1.0 


Carina 


20.10 ± 0.04 


0.025 


-1.73 ± 0.35 


-1.73 


Fornax 


20.74 ± 0.07 


0.025 


-0.99 ± 0.44 


-0.3 to -0.8 " 


Leo I 


22.07 ± 0.07 


0.03 


-1.35 ± 0.24 


-1.35 



4.2 explains these revisions in Sculptor and Fornax. 



TABLE 2 
Distances to the Galaxies 



Galaxy 


Distance 
modulus 


Ref. 


Sculptor 


19.71 ± 0.10 
19.64 ± 0.04 
19.67 ± 0.12 


Kaluzny et al. (1995) 
Rizzi et al. (2007a) 

Pictryzynski ct al. (2008) 


Garina 


20.09 ± 0.06 
20.06 ± 0.12 
20.19 ± 0.12 
20.11 ± 0.13 


Smecker-Hanc ct al. (1994) 

Mateo et al. (1998) 

DairOra et al. (2003) 

Pietryzynski et al. (2009) (avg. of J, K) 



Fornax 20.76 ± 0.10 Buoiiamio ct al. (1999) 

20.70 ± 0.12 Savianc ct al. (2000) (tip of RGB) 

20.76 ± 0.04 Savianc ct al. (2000) (HE) 

20.65 ±0.11 Bersier (2000) 

20.86 ± (0.04) Pietryzynski et al. (2003) 

20.66 ± (0.04) Mackey & Gilmore (2003) 
20.64 ± 0.09 Greco et al. (2007) 
20.72 ± 0.04 Rizzi et al. (2007b) 

20.75 ± 0.19 GuUieuszik et al. (2007) (tip of RGB) 

20.75 ± 0.11 GuUieuszik et al. (2007) (red clump) 

20.84 ± 0.14 Pietryzynski et al. (2009) 

Leo I 22.18 ± 0.11 Lee et al. (1993) 

22.00 ± 0.15 Caputo et al. (1999) 

22.04 ± 0.14 Held ct al. (2001) 

22.05 ± 0.18 Mcndcz ct al. (2002) 
22.02 ± 0.13 Bellazzini et al. (2004) 
22.04 ± 0.11 Held et al. (2010) 

paper; all mctallicities are on the CG97 scale froni Car- 
retta & Gratton 1997). Similarly, Kirby et al. (2009) find 
<[Fe/H]> = -1.58 ± 0.41 for a large sample within 10' 
of the center, while Helmi et al. (2006) find <[Fc/H]> = 
— 1.82 ± 0.34 for a sample out to the tidal radius (76f5; 
our <[Fe/H]> calculation based on their data). 

Revaz et al. (2009) explain this dichotomy with a 
young, metal-rich population (age < 2 Gyr), and an old 
metal-poor population (age > 9 Gyr). Their models in- 
dicate no intermediate-age population. Our carbon stars 
are more likely to belong to the younger population (see 
§ 4.2), which corresponds to the population dominating 
the core of the galaxy. To approximate the metallicity 
in the core, we have averaged the inner sample defined 
by Tolstoy et al. (2004) and the sample of Kirby et al. 
(2009), weighting by sample size (97 and 393, respec- 
tively), to arrive at an estimated <[Fe/H]> of —1.56 ± 
0.40. However, this value is too metal-poor compared to 
the metallicity of the younger population according to 
the models of Revaz et al. (2009), leading us to revise 
the metallicity in § 4.2. 

2.2.2. Carina 

Carina has experienced multiple, discrete star- 
formation events (Mighell 1990). Smecker-Hane et al. 



(1996) detected three main-sequence turn-off points, 
along with three locations for helium-burning stars (two 
HBs and a red clump projected onto the RGB), which 
Hurley-Keller et al. (1998) dated to three star-formation 
events early in the galaxy's history and ^7 and ^3 Gyr 
ago. This basic scenario has stood the test of another 
decade of observations (e.g., Monelli et al. 2003; Bono et 
al. 2010, and references therein). Models by Revaz et al. 
(2009) suggest that the intermediate population formed 
in a series of several bursts. 

Despite this complex star-formation history. Carina 
shows a relatively shallow metallicity gradient (Walker 
et al. 2009) and a MDF no broader than that of Sculp- 
tor. CaT observations with VLT/FLAMES by Helmi et 
al. (2006) cover Carina out to ~36', and from their data 
we determine that <[Fe/H]> = -1.81 ± 0.31. Other re- 
cent estimates are -1.72 ± 0.39 (Koch et al. 2006), -1.69 
± 0.51 (Koch et al. 2008), and -1.70 ± 0.19 (Bono et 
al. 2010). The mean of these results is -1.73 ± 0.35 (we 
have averaged the available uncertainties). 

2.2.3. Fornax 

Fornax has experienced a steadier rate of star forma- 
tion over the past several Gyr compared to Sculptor and 
Carina, resulting in a broader MDF. The CMD of For- 
nax has a wide RGB (Demers et al. 1979), which requires 
a range of mctallicities and/or ages. Fornax contains 
many carbon stars, proof of a substantial population of 
intermediate-age stars (Demers & Kunkel 1979; Aaron- 
son & Mould 1980). Further study made it apparent that 
episodes of star formation have continued to within the 
last few hundred Myr (Aaronson & Mould 1985; Buo- 
nanno et al. 1985, 1999), and that younger stars are 
more concentrated in the core of the galaxy (Stetson et 
al. 1998). More recent work has largely confirmed these 
earlier findings (e.g., Tolstoy et al. 2001; Pont et al. 2004; 
Helmi et al. 2006). 

CaT observations in large samples of RGB stars by 
Battaglia et al. (2006) trace the metallicity gradient in 
Fornax. They find a MDF within 24' of the center 
with <[Fe/H]> = -0.99 ± 0.44, compared to -1.52 ± 
0.46 outside 42' of the center (the quoted quantities are 
our determinations from their data). A study of CMDs 
within Fornax by Coleman & de Jong (2008) support 
the spectroscopic results. Our sample of carbon stars 
mostly conforms to the innermost sample considered by 
Battaglia et al. (2006), and we will adopt —0.99 as a 
starting metallicity for consideration. 

However, models by Revaz et al. (2009) suggest that 
the carbon stars could be significantly more metal-rich; 
most of the stars formed in the past few Gyr should have 
[Fe/H] in the range from ~ —0.3 to —0.8, which would 
make these carbon stars more similar to those in the 
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LMC and SMC than in the other three dwarf spheroidal 
galaxies in our sample. We return to this point in § 4.2 
below. 

2.2.4. Leo I 

Leo I is the most distant of the dwarf galaxies in the 
Milky Way system of the Local Group. In fact, it is 
unclear whether or not it is gravitationally bound to the 
Local Group (e.g., Lepine et al. 2011). Its distance and 
its proximity to Regulus have made observations more 
challenging than for the other dwarfs considered here. 
Nonetheless, a picture has emerged of a galaxy with the 
contradictory properties of a relatively young population 
and relatively metal-poor abundances (Lee et al. 1993). 
The majority of the visible stars in the galaxy appear 
to have formed 3-7 Gyr ago (Demers et al. 1994) or 1-7 
Gyr ago (Gallart et al. 1999). 

Despite the ongoing star formation, gradients in the 
population are subtle (GuUieuszik et al. 2009; Held et al. 
2010), and the metallicity shows a fairly narrow and well- 
defined distribution. CaT spectra give <[Fe/H]> values 
of -1.34 ± 0.26 (102 stars; Bosler et al. 2007), -1.31 ± 
0.25 (58 stars; Koch et al. 2007), and -1.41 ± 0.21 (54 
stars; GuUieuszik et al. 2009). Combining these results 
and weighting by their sample size yields <[Fe/H]> = 
-1.35 ± 0.24. 

2.3. Stellar samples 

The carbon stars in our sample were observed in two 
Spitzer programs. The first, a Cycle 2 program, included 
five carbon stars in Fornax (published by Matsuura et al. 
2007) and three carbon-star candidates in Leo I. The sec- 
ond program followed in Cycle 3 and included six carbon 
stars in Fornax, three in Carina, and two in Sculptor. 
Sloan et al. (2009) published one of the two Sculptor 
spectra. Table 3 gives the names, positions, and NIR 
fluxes from 2MASS (Skrutskie et al. 2006) of the stars in 
our sample. Figure 1 shows where our targets are located 
in each galaxy. 

Fornax is well known as an abundant source of carbon 
stars, starting with the initial detection of several candi- 
dates by Demers & Kunkel (1979) and the spectroscopic 
confirmation of six by Aaronson & Mould (1980). By 
1999, the number of confirmed carbon stars had climbed 
to 104 (Azzopardi et al. 1999). Multiple programs have 
searched Fornax for LPVs. Demers & Irwin (1987) found 
30 candidates, but no Mira variables. Bersier & Wood 
(2002) identified 85 candidates, but did not attempt to 
determine periods. Whitelock et al. (2009) published the 
results of a thorough NIR monitoring program from the 
SAAO. The Cycle 2 program selected five targets in For- 
nax based on the SAAO observations, although at that 
time the mean NIR magnitudes and periods were not 
known. The six targets in Cycle 3 were selected based 
on 1-5 /im spectra and photometry with the Infrared 
Spectrometer and Array Camera (ISAAC) on the VLT. 
Groenewegen et al. (2009a) pubhshed the results of these 
observations, which include eight of our 11 targets. The 
combination of the various observations confirms that all 
of our sample are carbon-rich. The SAAO photometry 
provides periods for six of our 11 targets in Fornax. Four 
of these are Mira variables, while two are semi-regulars. 

Our search for targets in Sculptor began with the 
2MASS survey. Two targets fulfilled our criteria that 
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Fig. 1. — The locations of our targets within each galaxy. The 
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TABLE 3 

Targets 



Source* 


Position (J2000) 




2MASS photometry 




Other 


name 


RA Dec. 


J 


H 


Ks 


designations* 



MAG 29 


00 


59 


53. 


.67 


-33 


38 


30. 


.8 


14. 


.846 


± 


0. 


.038 


13. 


.144 


± 


0.031 


11. 


.603 


± 


0.021 




Scl V78 V544 


00 


59 


58. 


.94 


-33 


28 


35. 


.2 


13. 


.399 


± 


0. 


.021 


12. 


.633 


± 


0.025 


12. 


.273 


± 


0.023 


ALW Scl 3 


For BW 2 


02 


38 


06. 


,19 


-34 


31 


19, 


.4 


16, 


.052 


± 


0. 


,104 


14, 


.483 


± 


0.055 


13, 


.315 




0.048 


GLM 31 


For BTH 13-23 


02 


38 


50. 


,56 


-34 


40 


32, 


.0 


16, 


.106 


± 


0. 


,094 


14. 


.525 


± 


0.053 


12, 


.879 


± 


0.029 




For BTH 12-4 


02 


39 


12. 


,33 


-34 


32 


45, 


.0 


14, 


.722 


± 


0. 


,033 


13. 


.262 


± 


0.038 


12, 


.120 


± 


0.024 


GLM 25 


For BTH 3-129 


02 


39 


41. 


,60 


-34 


35 


56, 


.7 












15. 


.970 


± 


0.205 


14, 


.164 


± 


0.070 




For DK 18 


02 


39 


54. 


,21 


-34 


38 


36, 


.9 


15, 


.601 


± 


0. 


,064 


14, 


.162 


± 


0.049 


13, 


.167 


± 


0.033 


GLM 24 


For DK 52 


02 


40 


06. 


,66 


-34 


23 


22. 


.3 


14. 


.485 


± 


0. 


,029 


13. 


.377 


± 


0.034 


12. 


.618 


± 


0.027 


DDB 17, GLM 13 


For DI 2 


02 


40 


09. 


,47 


-34 


06 


25. 


.7 


15. 


.790 


± 


0. 


,075 


14, 


.556 


± 


0.068 


13. 


.668 


± 


0.052 


GLM 16 


For WEL CIO 


02 


40 


10. 


,17 


-34 


33 


21. 


.9 


14. 


.063 


± 


0. 


,026 


13. 


.122 


± 


0.021 


12. 


.545 


± 


0.029 


DI 20, SHS 105, BW 62, DDB 


For BW 69 


02 


40 


17. 


,79 


-34 


27 


35. 


.8 


15. 


.424 


± 


0. 


,063 


14. 


.122 


± 


0.049 


13. 


.182 


± 


0.035 


GLM 21 


For BW 75 


02 


40 


31. 


,23 


-34 


28 


44. 


.2 


14. 


.745 


± 


0. 


,043 


13. 


.689 


± 


0.046 


13. 


.072 


± 


0.037 


DDB 22, BTH 6-13, GLM 17 


For BW 83 


02 


41 


03. 


,56 


-34 


48 


05, 


.4 


14, 


.441 


± 


0. 


,035 


13, 


.365 


± 


0.034 


12, 


.694 


± 


0.034 


DDB 25, MAG 30, GLM 27 


ALW Car 2 


06 


41 


13. 


,53 


-50 


54 


25, 


.0 


13, 


.925 


± 


0. 


,024 


13. 


.073 


± 


0.028 


12, 


.658 




0.027 




Car MCA C3 


06 


41 


41. 


,45 


-50 


58 


08, 


.1 


13, 


.742 


± 


0. 


,022 


12. 


.805 


± 


0.023 


12, 


.340 


± 


0.026 


ALW Car 6 


Car MCA C5 


06 


42 


10. 


,35 


-50 


56 


24. 


.0 


13. 


.940 


± 


0. 


,023 


13. 


.159 


± 


0.027 


12. 


.785 


± 


0.029 


ALW Car 10 



Leo I MFT C 10 08 22.25 +12 17 57.1 ■•• 16.195 ± 0.214 14.225 ± 0.060 HGR 8717 

Leo I MFT A 10 08 29.28 +12 18 51.6 17.134 + 0.202 15.429 + 0.103 14.025 + 0.053 HGR 6343 
Leo I MFT E 10 09 00.5 +12 19 01 

=1 V78 = van Agt (1978); DK = Dcmcrs & Kunkcl (1979); MCA = Mould ct al. (1982); ALW = Azzopardi ct al. (1985. 1986); WEL = Wcstcrlund 
et al. (1987); DI = Dcmcrs & Irwin (1987); SHS = Stetson ct al. (1998); BW = Bcrsicr & Wood (2002); DDB = Dcmcrs ct al. (2002); MFT = 
Menzies et al. (2002); MAG = Mauron et al. (2004); BTH = Battaglia et al. (2006); GLM = Groenewegen et al. (2009a); HGR = Held et al. 
(2010). 



TABLE 4 
Optical photometry 



Source 




Photometry* 








Ref.'^ 


name 


B 


V 




R 




I 




MAG 29 






20.22 






18.04 


USNO-B 


Scl V78 V544 


20.20 


16.55 


16.53 


+ 


0.49 


16.70 


USNO-B 


For BW 2 




20.23 + 0.07 


19.21 


+ 


1.58 


16.47 + 0.05 


BW, USNO-B 


For BTH 13-23 














For BTH 12-4 
















For BTH 3-129 
















For DK 18 


21.45 


18.57 + 0.60 


17.51 


+ 


1.46 


16.01 


DK, USNO-B 


For DK 52 


22.06 


19.74 + 0.33 


17.93 


+ 


0.51 


18.15 


DK, WEL, USNO-B 


For DI 2 


21.4 


18.3 + 0.7 


17.07 


+ 


0.42 


16.75 


WEL, DI, USNO-B 


For WEL CIO 


22.31 + 0.03 


19.41 + 0.05 


18.24 


+ 


1.08 


15.73 + 0.04 


SHS, BW, USNO-B 


For BW 69 


21.93 


19.99 + 0.05 


19.80 


+ 


0.74 


16.64 + 0.04 


BW, USNO-B, GSC 


For BW 75 


20.63 


20.06 + 0.03 


18.14 


+ 


0.33 


16.99 + 0.07 


GSC. BW. USNO-B 


For BW 83 


23.26 


20.41 + 0.10 


17.63 


+ 


0.98 


16.62 + 0.06 


BW, USNO-B, NOMAD 


ALW Car 2 


18.89 ± 0.64 


17.64 


16.37 


± 


0.01 


15.64 


USNO-B, NOMAD, K06 


Car MCA C3 


18.43 


17.48 


15.86 






15.27 


USNO-B, NOMAD 


Car MCA C5 


18.95 + 0.71 


17.26 + 0.09 


16.01 


+ 


0.33 


15.11 


M82, USNO-B, NOMAD 



Leo I MFT C • • • 
Leo I MFT A • • • 
Leo I MFT E • • • 

* Entries in bold are mean magnitudes and their standard deviation or amplitude from multiple photometric 
observations. If only the mean magnitude is bold, then it is followed by an uncertainty in the mean. 

DK = Dcmcrs & Kunkcl (1979); M82 = Mould ct al. (1982); WEL = Wcstcrlund ct al. (1987); DI = Dcmcrs 
& Irwin (1987); SHS = Stetson et al. (1998); BW = Bersier & Wood (2002); USNO-B = Monet et al. (2003); 
NOMAD = Zacharias et al. (2004); K06 = Koch et al. (2006); GSC = Lasker et al. (2008). 
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J - Ks > 1.1 and Ks < 12.5. V78 V544 was originally 
identified as a LPV by van Agt (1978). This star had 
the reddest J — K color (1.19) in the study by Frogel et 
al. (1982), but they were unable to determine if it was 
carbon-rich spectroscopically. Azzopardi et al. (1986) 
made that confirmation, noting that their list of eight 
carbon stars was likely to be complete. However, the 
optical surveys in use at that time missed the sources 
embedded in the most optically thick shells. Our other 
Sculptor target, MAG 29, has J - = 3.24. Mauron 
et al. (2004) first noticed this source in their search of 
2MASS targets in the direction of several dwarf galax- 
ies, but they estimated its distance to be ^50 kpc, in 
the foreground of the Sculptor dwarf. Sloan et al. (2009) 
published an early version of the IRS spectrum of MAG 
29. Using two different infrared color-magnitude rela- 
tions, they estimated its distance to be 84 zb 13 kpc, 
consistent with the distance of Sculptor in Table 1, 85 ± 
2 kpc. Ncar-infrarcd spectra of V78 V544 and MAG 29 
by Grocnewcgen et al. (2009a) confirm their carbon-rich 
nature. Menzics ct al. (2011) find that both are Mira 
variables and that the P-L relation gives distances for 
these two consistent with membership in Sculptor. 

We also used the 2MASS survey to search Carina 
for suitable targets. Four sources fulfilled the crite- 
ria J — Kg > 1.1 and Kg < 13.0 (excluding obvious 
foreground sources), and of these we observed three^^. 
Mould et al. (1982) originally identified two of our three 
targets as carbon stars, and wc adopt their names for 
them. Azzopardi ct al. (1986) included all four of the red 
sources in their list of nine spectroscopically confirmed 
carbon stars in Carina. They stated that the list should 
be complete for the area observed. 

More carbon stars have been detected in Leo I than 
in Sculptor or Carina. Azzopardi et al. (1986) listed 16 
spectroscopically confirmed carbon stars and two candi- 
date carbon stars in Leo I. However, none of these are 
very red. A more recent NIR survey from the SAAO de- 
tected five highly reddened stars with J — Kg > 2 (Men- 
zies et al. 2002). In our Cycle 2 program, we observed 
the three reddest, all with J — Kg > 3. Mcnzies et al. 
(2010) determined periods for our three Leo I targets as 
part of a larger efli'ort which identified 26 AGB variables 
in the galaxy from the SAAO. While it is quite likely 
that all three of our targets are carbon-rich, the spec- 
tra presented here are our first chance to confirm their 
chemistry. 

Table 4 gives the best available optical photometry for 
the sources in our sample. Entries in bold are average 
magnitudes and the standard deviations when the quoted 
sources give multiple measurements or provide a mean. 
Table 5 presents variability classes, periods, and NIR 
photometry. The entries in bold in this table are mean 
magnitudes and peak-to-peak amplitudes published by 
the SAAO (references are given in the table notes). 

3. OBSERVATIONS AND DATA REDUCTION 

The IRS obtained infrared spectra of our 19 targets in 
the standard staring mode, observing each target in the 
two Short-Low (SL) apertures, producing spectra from 5 
to 14 /xm. SL order 2 (SL2) covers the 5.1-7.5 /xm re- 
gion, while SL order 1 (SLl) covers 7.5-14.2 /xm. Two 

The unobserved source is ALW Ccir 7. 



relatively bright targets were also observed in the two 
Long-Low (LL) apertures (LL2: 14-20.5 ^m; LLl: 20.5- 
37 ^m). Table 6 presents the details for each observa- 
tion. The SL2 and LL2 spectra include a short piece of 
a first-order spectrum, the "bonus" order. These bonus 
orders provided overlap between SL2 and SLl (and LL2 
and LLl), making it possible to determine multiplicative 
corrections to remove discontinuities between segments 
which arise from pointing shifts during the IRS integra- 
tions. Each source was observed in two separate nod 
positions in each aperture, requiring four separate point- 
ings for the SL-only spectra and eight for the two which 
included LL. 

Our data analysis began with the fiatfickled imagcs^^ 
produced by the S18.18 pipeline at the Spitzer Science 
Center (SSC). Before extracting spectra from images, we 
removed the background by subtracting the correspond- 
ing image with the source in a different position, either 
in the other nod position in the same aperture (a nod 
difference) or in the same nod position, but in the other 
aperture (an aperture difference). The SL images in Cy- 
cle 3 used aperture differences, but the Cycle 2 images 
required nod differences, due to the mismatching num- 
ber of observations in SL2 and SLl. Nod differences were 
used for the two LL spectra. 

The differenced IRS images were then cleaned, using 
the IMCLEAN procedure, which is similar to the irsclean 
procedure available from the SSC. Pixels were replaced 
with an average computed from surrounding rows if they 
were flagged as bad or if they were included in the cam- 
paign masks of rogue pixels distributed by the SSC. We 
generally treated a pixel as a rogue if it had been flagged 
as such twice in the current or any prior IRS campaign. 
The number of rogue pixels steadily grew over the course 
of the cryogenic Spitzer mission, and not all rogue pix- 
els were ever flagged. Wc added several additional pix- 
els to the rogue masks for the Cycle 3 data when we 
could see the impact of consistently misbehaving pix- 
els on our final spectra. This step was crucial for im- 
proving the signal/noise (S/N) ratio of our faint spec- 
tra, because these unflaggcd rogue pixels contribute non- 
gaussian noise which becomes more significant at low flux 
levels. 

Our extraction of spectra from the images relied on the 
optimal extraction algorithm developed at Cornell and 
described in detail by Lebouteiller et al. (2010). This al- 
gorithm fits a wavelength-dependent point-spread func- 
tion (PSF) to each row of the spectral image, reducing 
the impact of noise from pixels containing little fiux from 
the source. For point sources, it substantially improves 
the S/N compared to spectra extracted with more con- 
ventional algorithms. 

Spectra from individual images with the source in a 
given nod position were then co-added. Wc used a spike- 
rejection algorithm to fiirther reduce the effect of non- 
gaussian noise components when combining the spectra 
from the two nod positions. At that stage, we re-assessed 
the propagated noise. If the uncertainty as measured by 
comparing the two nod positions was larger, we used this 
value instead. Finally, spectra from different spectral or- 
ders were combined using a "stitch-and-trim" algorithm, 
first applying multiplicative corrections to remove dis- 

Basic calibrated data, or BCD files. 
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TABLE 5 
Variability 



Source 


Var. 


Period 








Photometry^ 








Ref.'' 


name 


class 


(days) 


J 


AJ 


ll 




Ks 






A /T A on 

MACj z9 


Mira 


554 


14.35 






12.94 




11.44 





87 


Mil 


ocl V7» V544 


Mira 


189 


13.78 






12.90 




12.38 





42 


Mil 


For BW 2 


var. 




16.05 ± 0.10 






14.48 ± 0.06 




13.32 ± 0.05 






BW, 2MASS 


^or BiH 13-23 


Mira 


350 


17.09 


1 


28 


15.33 


1.24 


13.63 


1 


02 


WOO 


I'or BiH iz-4 


Mira 


470 


16.01 


1 


44 


14.31 


1.14 


12.91 





96 


W09 


For BTH 3-129 


Mira 


400 


18.00 


1 


64 


15.83 


1.43 


13.90 


1 


16 


W09 


For DK 18 


vax. 




14.71 


1 


00 


13.59 


0.76 


12.78 





43 


WOO 


tor DK 52 


var. 




14.76 





61 


13.60 


0.52 


12.80 





32 


WOO 


For DI 2 


irr. 




15.79 ± 0.07 






14.56 ± 0.07 




13.67 ± 0.05 






Dl, 2MAbt3 


For WEL CIO 


SR 


317 


15.09 





64 


13.87 


0.43 


13.05 





26 


DI, WOO 


For BW 69 


SR 


340 


16.24 






14.85 




13.61 






WOO 


For BW 75 


var. 




15.12 


1 


00 


13.94 


0.75 


13.16 





51 


WOO 


For BW 83 


Mira 


280 


14.87 





69 


13.77 


0.56 


12.99 





52 


W09 


ALW Car 2 


var. 




13.93 ± 0.03 






13.07 ± 0.03 




12.66 ± 0.03 






SAAO, 2MASS<= 


Car MCA C3 


var. 




13.74 ± 0.03 






12.81 ± 0.02 




12.34 ± 0.03 






SAAO, 2MASS<= 


Car MCA C5 


var. 




13.94 ± 0.03 






13.16 ± 0.03 




12.79 ± 0.03 






SAAO, 2MASS<= 


Leo I MFT C 


Mira 


523 


17.46 


1 


52 


15.68 


1.29 


13.98 


1 


03 


MIO 


Leo I MFT A 


Mira 


336 


17.62 


1 


23 


15.88 


1.01 


14.39 





81 


MIO 


Leo I MFT E 


Mira 


283 


19.18 


1 


87 


17.25 


1.23 


15.63 


1 


17 


MIO 



^ Entries in bold are mean magnitudes and pcak-to-pcak amplitudes based on light-curve analysis. 

DI = Demers & Irwin (1987); BW = Bersier & Wood (2002); 2MASS = Skrutskie et al. (2006); WOO = Whitelock et al. 
(2009); MIO = Menzies et al. (2010); Mil = Mcnzies ct al. (2011); SAAO = Unpublished communication from SAAO. 

The SAAO identifies the target as variable; the photometry are from 2MASS. 



TABLE 6 

IRS OnSER\"ATIOKS 



Source 


AOR 


Program 




Observing date 




On-source ir 


name 


key 




day 


JD-2400000.5 


SL2 


SLl 


MAG 29 


18050816 


30333 


2006 


Dec 19 


54088 


9 


240 


240 


Scl V78 V544 


18051072 


30333 


2006 


Dec 20 


54089 


1 


2880 


5280 


For BW 2 


18053376 


30333 


2007 


Feb 08 


54139 


7 


480 


480 


For BTH 13-23 


14540544 


20357 


2006 


Jan 30 


53765 


6 


2880 


3360 


For BTH 12-4 


14541056 


20357 


2006 


Jan 27 


53762 


9 


2880 


3360 


For BTH 3-129 


14540800 


20357 


2006 


Jan 30 


53765 


5 


2880 


3360 


For DK 18 


18052864 


30333 


2007 


Feb 08 


54139 


6 


1440 


1440 


For DK 52 


18052096 


30333 


2006 


Dec 10 


54088 


9 


1440 


1440 


For DI 2 


18052352 


30333 


2007 


Feb 06 


54137 


6 


2880 


2880 


For WEL CIO 


14541824 


20357 


2006 


Jan 27 


53762 


9 


2880 


3360 


For BW 69 


18052608 


30333 


2007 


Feb 07 


54138 


9 


1440 


1440 


For BW 75 


14542080 


20357 


2006 


Jan 27 


53762 


8 


2880 


3360 


For BW 83 


18053120 


30333 


2007 


Feb 08 


54139 


7 


1440 


1440 


ALW Car 2 


18051584 


30333 


2007 


Mar 12 


54171 


9 


2880 


5280 


Car MCA C3 


18051328 


30333 


2007 


Mar 12 


54171 


8 


2880 


4800 


Car MCA C5 


18051840 


30333 


2007 


Mar 13 


54172 





3840 


5760 


Leo I MFT C 


14545152 


20357 


2006 


May 25 


53880 





2880 


3360 


Leo I MFT A 


14545408 


20357 


2006 


May 25 


53880 


1 


2880 


3360 


Leo I MFT E 


14545920 


20357 


2006 


May 25 


53880 


2 


2880 


3120 



LL2 



LLl 



960 



2880 



960 



2880 



continuities between spectral segments, then truncating 
invalid data from the ends of each segment. The correc- 
tions were typically on the order of -^5%, although they 
could be as large as 15%. 

The photometric calibration of the spectra has changed 
slightly from previous publications from the IRS team at 
Cornell, as outlined by Lcbouteillcr ct al. (2011). We 
use HR 6348 (KO III) as the standard for SL, as before, 
but the assumed truth spectrum for this source has been 
shifted down 5% to align with the updated calibration of 
the Multiband Imaging Photometer for Spitzer (MIPS) 



at 24 ^m (Riekc et al. 2008). We used HR 6348 and 
HD 173511 (K5 III) for LL, with the latter spectrum 
shifted similarly. 

Figures 2-4 present the final SL spectra for our 19 tar- 
gets. Figure 5 presents the full spectra for the two sources 
also observed in LL. 

4. REFINING THE METALLICITY ESTIMATES 

4.1. Bolom,etric m,agnitudes 

We estimate bolometric magnitudes by integrating the 
IRS spectra, combined with the available optical and 
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TABLE 7 

BOLOMETRIC MAGNITUDES 



Source 


Our 


External 


Other values for Mur,i° 


name 


Mboi 


uncertainty'' 


M07 


L08 


G09 


SAAO 




— 4.9U ± (J.U4 












bcl V ( o Vo44 


—4.20 ± U.U4 












For BW 2 


-4.41 ± 0.29 


0.25 






-4.05 




ror r>iH lo-zo 


—4.28 ± U.U4 


0.30 


-4.90 


-4.84 




-4.45 


tor D iri 12-4 


— 4.yU ± U.U4 


0.29 


—5.52 


-5.26 


-5.20 


—4.80 


For BTH 3-129 


-4.65 ± 0.04 


0.17 


-4.95 






-4.68 


For DK 18 


-4.88 ± 0.04 


0.54 






-4.11 




tor UK 52 


— 4.62 ± 0.04 


0.02 






-4.65 




tor UL 2 


—4.06 ± 0.29 


0.34 






-3.58 




For WEL CIO 


-4.61 ± 0.04 


0.20 


-4.60: 






-4.21 


For BW 69 


-3.98 ± 0.04 


0.06 






-4.07 




For BW 75 


-4.35 ± 0.04 


0.22 


-4.67 




-4.25 




For BW 83 


-4.57 ± 0.04 


0.19 






-4.61 


-4.27 


ALW Car 2 


-4.54 ± 0.29 












Car MCA C3 


-4.82 ± 0.29 












Car MCA C5 


-4.57 ± 0.29 












Leo I MET C 


-4.77 ± 0.13 


0.47 








-5.44 


Leo I MET A 


-4.43 ± 0.13 


0.23 








-4.75 


Leo I MET E 


-4.34 ± 0.13 


0.42 








-3.74 



^ The standard deviation of all of the given values for Mboi . 

^ Adjusted to our adopted distance modulus. M07 = Matsuura et al. (2007); LOS 
= Lagadec et al. (2008); G09 = Groenewegen et al. {2009a), SAAO = Whitelock et 
al. (2009) for Fornax and Menzies et al. (2010) for Leo I. 



For BTH 13-23- 6.5 mjy 




Fig. 2. — Optimally extracted SL spectra of the five Fornax car- 
bon stars originally presented by Matsuura et al. (2007). The fluxes 
measured at 6.0—6.5 fim are given beside each target name. 



NIR photometry. Table 7 presents the resuhs and com- 
pares them with previously published estimates. Over 
the region covered by the spectrum, simple integration 
suffices. Below 5.1 /im, we integrate on a grid of F\ vs. 
wavelength linearly interpolated through the photometry 
in Tables 4 and 5. We extrapolate with a Rayleigh- Jeans 
tail at the long-wavelength end and a Wien distribution 
at the short-wavelength end. The final magnitudes are 
scaled to the distances given in Table 1. 

This paper marks a shift from previous determinations 
of bolometric magnitude. Before, we did not consider 



For BW 2 - 5.1 mJy 




Fig. 3. — Optimally extracted SL spectra of the six new carbon 
stars in Fornax. Fluxes measured at 6.0—6.5 fim follow each target 
name. 



BVRI photometry. The change is significant. The Wien 
distribution extrapolated from J overestimates the flux 
in the visual regime by ~0.2 magnitudes for our bluest 
sources, because the combination of molecular band ab- 
sorption and dust extinction drops more quickly with de- 
creasing wavelength. The largest difference in our sam- 
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Leo 1 MFT E - 2.2 mJy 

_i I I I I I I I I I I I I L 



10 
A (miti) 



12 



14 



Fig. 4. — Optimally extracted SL spectra of the sources in Sculp- 
tor, Carina, and Leo I. The fluxes after each target name are mea- 
sured at 6.0-6.5 /im. 



jA MAG 29 - 


?\,y. For BW 2 

i ^>i,j 

'A/tiidijiAj L 111 











10 



15 



20 

A (mid) 



25 



30 



35 



Fig. 5.— The full IRS spectra for MAG 29 and For BW 2, which 
were observed in both SL and LL. The dotted lines mark the zero- 
flux levels for both spectra. 

pie is 0.27 magnitudes'^. The difference decreases to 
zero past J — A' ~ 3. A line can be fitted to this shift: 
AMboi = 0.27 - 0.075(J - K), although the scatter is 
^0.1 magnitudes. It would be appropriate to apply this 
correction to our previously published bolometric magni- 
tudes for carbon stars (e.g., Sloan et al. 2006, and most 
of the other IRS papers cited in §1). 

We have simplified the treatment of the BVRI pho- 
tometry by not considering differences in photometric 

15 The star was For DK 52. 



systems or attempting transformations among them. 
Such an effort would have virtually no impact on the 
bolometric magnitudes reported here. However, for the 
NIR photometry, we did distinguish between the SAAO 
and 2MASS systems, as it is in this wavelength regime 
that the spectral energy distributions (SEDs) of our 
sources peak'^. 

The impact of the correction for interstellar reddening 
is much smaller. It brightens the bluest stars by only 
^^0.02 magnitudes in our samples, with the correction 
decreasing to nearly zero for the most enshrouded carbon 
stars. The redder the J — K color, the more the peak of 
the SED has shifted away from wavelengths most affected 
by interstellar reddening. 

The most significant source of uncertainty in our sam- 
ple is the variability of the star. The JHK photometry 
dominates the final bolometric magnitude. Comparing 
the mean magnitudes from the SAAO to the 2MASS 
data reveals an average difference (|AMboi|) of 0.29 mag- 
nitudes, in either direction, with differences as large as 
0.6 magnitudes. Interestingly, we get similarly large dif- 
ferences if we only consider the sources identified as semi- 
regulars, irregulars, or simply "variables". For any vari- 
able star, the mean magnitude is definitely preferable to 
NIR photometry that only sparsely covers the period of 
variability. 

For the sources with mean NIR magnitudes, the un- 
certainties in bolometric magnitude in Table 7 are just 
the uncertainty in distance modulus. For the remaining 
five sources, we set the uncertainties to 0.29 to reflect 
the limitations of the 2MASS photometry. The reader 
should bear in mind, though, that errors as large as 0.6 
magnitudes are possible. 

Table 7 also lists what we describe as "external" un- 
certainties in A/boi- These are the standard deviation of 
all values of Mboi for a given star in Table 7. The values 
from other authors have been adjusted to our adopted 
distance moduli (in Table 1). The values published by 
Matsuura et al. (2007) are based on earlier versions of the 
IRS data for the five stars in common between this work 
and theirs. They used these data along with the NIR 
photometry from the SAAO and fitted radiative trans- 
fer models to determine the luminosity. Four of their 
bolometric magnitudes appear in their Table 4; we re- 
constructed the fifth from the luminosity given in the 
text. Lagadec et al. (2008) used 2MASS photometry 
and applied the bolometric corrections defined by White- 
lock et al. (2006), which were calibrated from Galactic 
carbon stars with photometry from the optical into the 
mid-infrared, including mean magnitudes at JHK, IRAS 
data to 25 fim and MS'X data to 15 /im'^. Groenewegen 
et al. (2009a) applied the bolometric corrections defined 
by BesseU & Wood (1984) to 2MASS photometry. The 
bolometric magnitudes from the SAAO in Table 7 are 
based on mean magnitudes at JHK and the bolometric 
corrections of Whitelock et al. (2006). The external un- 
certainties may well overestimate our actual uncertainty 
in A/boi as they reflect systematic differences between 

1^ Cohen, et al. (2003) deflned the the central wavelengths and 
zero-magnitude fluxes for 2MASS. For SAAO, Nagayama et al. 
(2003) defined the wavelengths, and we scaled the zero-magnitude 
fiuxes from the 2MASS data. 

^"^ IRAS: the Infrared Astronomical Satellite (Beichman et al. 
1988); MSX: the Mid-course Space Experiment (Egan et al. 2003). 
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Period-luminosity reiation - Wiiiteioci( et ai. (2009) ~ 

_j i I 1 I 1 I 1 I I I I 1 I lJ 

200 400 600 800 

Period (days) 

Fig. 6. — The 11 sources for which we have periods plotted on 
a period-luminosity (P-L) diagram and compared to evolutionary 
tracks by Vassiliadis & Wood (1993) and the P-L relation of White- 
lock et al. (2009). The masses given are initial values in solar 
masses. For Mboi, the error bars are the larger of the uncertainties 
in Table 7, with a minimum assumed uncertainty of 0.15 magni- 
tudes (which affects Fornax BW 69 and the two Sculptor objects). 
We have assumed a 10% uncertainty in period. To avoid clutter 
in the figure, we have shortened the source names. The two open 
circles mark the semi-regular variables. 



TABLE 8 
Possible stellar ages 



Source 


Estimated 


Estimated 


name 


initial mass (Mq) 


age (Gyr) 


MAG 29 


1.5-1.9 


0.9-1.7 


Scl V78 V544 


1.2-1.8 


1.0-3.5 


For BTH 13-23 


0.97-1.5 


1.8-11 


For BTH 12-4 


1.4-2.2 


0.6-2.7 


For BTH 3-129 


1.3-1.7 


1.2-3.3 


For WEL CIO 


1.3-1.9 


0.9-3.3 


For BW 69 


0.94-1.1 


5.5-12 


For BW 83 


1.3-2.0 


0.8-3.3 


Leo I MET C 


1.0-2.3 


0.5-6.9 


Leo I MET A 


1.0-1.5 


1.7-6.9 


Leo I MET E 


0.98-1.9 


0.8-7.3 



groups with access to different data, but including these 
systematics will force us to be cautious with our Mboi 
results. 

4.2. Masses and metallicities 

We have pulsation periods for 11 of our 19 targets, in- 
cluding both carbon stars in Sculptor, all three in Leo I, 
and six of the 11 in Fornax (see Table 5). Figure 6 com- 
pares their periods and bolometric magnitudes to evolu- 
tionary tracks by Vassiliadis & Wood (1993), allowing us 
to make rough estimates of their initial mass, and thus 
their age. The evolutionary tracks convert the relatively 
robust period-mass-radius relation from stellar pulsation 
theory into a period-mass-luminosity relation by assum- 
ing a relation between effective temperature and luminos- 
ity. The mass-loss history must be estimated to convert 
current to initial mass, and the histories utilized by Vas- 
siliadis & Wood (1993) have successfully reproduced the 
luminosities at the tips of the RGB in Magellanic clus- 
ters. More theoretical approaches to the problem (e.g.. 



Kamath et al. 2011) basically validate the older evolu- 
tionary tracks. 

The error bars in Figure 6 are the larger of the internal 
and external uncertainties for Mboi in Table 7, although 
we have assumed a minimum uncertainty of 0.15 magni- 
tudes, which affects MAG 29 and V78 V544 in Sculptor 
and For BW 69. For period, we assumed a 10% uncer- 
tainty, based on the fact that Groenewegen et al. (2009b) 
have found shifts in period on this order for some of the 
sources in their sample (not counting those sources with 
possible changes in pulsation mode). Our conservative 
approach to the uncertainties leads to a fairly wide range 
of initial masses and ages in our sample, but even these 
broad limits can help constrain the metallicities. 

Table 8 gives our rough estimates of initial masses and 
current ages, based on Figure 6. Assuming that the stel- 
lar lifetimes scale as M~^-^ would overestimate the ages, 
because more metal-poor stars evolve off of the main se- 
quence more quickly. To estimate the ages, we have used 
models by Maraston (2005), who give main-sequence 
turn-off masses at metallicities of [Fe/H] = —1.35, —0.33, 
and 0.0. The first metallicity matches Leo I perfectly. 
For Sculptor we spline interpolated to [Fe/H] = —1.0. 
For Fornax, we spline interpolated to —0.3 and —0.8, us- 
ing the former for the lower bound on the mass and the 
latter for the upper bound. To get from main-sequence 
turn-off point to the AGB, we assumed that a star spends 
a time on the red giant branch equal to 10% of its main- 
sequence lifetime. 

The likely ages of both stars in Sculptor are consistent 
with the models of Revaz et al. (2009), which show vir- 
tually no star formation in Sculptor over the period from 
~8 to ~2 Gyr ago. Thus both stars must be younger 
than ~2 Gyr old. Their models indicate that stars pro- 
duced in this recent round of star formation have [Fe/H] 
~ —1.0. This value is significantly higher than the mean 
in Table 1 and the value of —1.4 assumed by Sloan et al. 
(2009). 

Menzies et al. (2011) also placed MAG 29 in the 
younger population in Sculptor, with an age of ^-^1-2 Gyr. 
However, they made this estimate based on its pulsa- 
tion period, which they argue is a good diagnostic for 
age (see their § 6 and references therein. They assigned 
Scl V78 V544 in the older population due to its shorter 
pulsation period, while we find that its luminosity is more 
consistent with the younger population (age < 2 Gyr). 
The discrepancy between the two approaches arises from 
the significant spread in periods possible for stars of a 
given mass. 

Two of the targets in Fornax, WEL CIO and BW 69, 
are semi-regular pulsators. The location of For BW 69 
amongst the other data suggests that it is a fundamental- 
mode pulsator. If it were pulsating in the first overtone, 
then its position in Figure 6 would correspond to a pe- 
riod ^2.2 times longer (Wood & Sebo 1996), which would 
imply an unreasonably low mass. If For WEL CIO were 
pulsating in the first overtone, it would shift to the right- 
most point in the figure and become something of an out- 
lier. While we suspect that it is pulsating in the funda- 
mental mode, an error here would have only a small im- 
pact on its estimated mass, because the track for 1.5 Mq 
is nearly horizontal. 

Four of the six sources with periods in Fornax have 
relatively young ages of ~3.3 Gyr or less. The models of 
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Revaz et al. (2009) suggest stars of this age would have 
corresponding metalhcities in the range —0.8 < [Fc/H] 
< —0.3. This range is fairly constant over the likely time 
frame, and it is more similar to the Magellanic samples 
than to the other dwarf spheroidals considered here. The 
uncertain mass of BTH 13-23 leads to an unconstrained 
age and a poorly constrained mctallicity. BW 69 looks to 
be at least 5.5 Gyr old. Using the models of Revaz et al. 
(2009) as a guide, its maximum mctallicity is ~— 0.5, the 
mean mctallicity of the SMC, although it cold be more 
metal-poor. 

Five stars in Fornax without periods do not appear in 
Figure 6. Two were outside the survey area of Whitelock 
ct al. (2009), and both were identified as variables by 
other observers (see Table 5). Whitelock et al. (2009) 
identified the other three as variables but were unable to 
report a period. All five bolometric magnitudes do little 
to constrain their likely masses, ages, and metalhcities, 
which we will assume are similar to the sources for which 
we have periods. 

The three carbon stars in Leo I are probably younger 
than ~7 Mq, but their relatively unconstrained masses 
allow us to say little more about their age. Fortunately 
for our efforts to constrain their mctallicity, the MDF for 
Leo I is narrow (<[Fc/H]> = -1.35 ± 0.24, sec § 2.2). 
No matter their mass and age, these stars arc likely to be 
more metal-poor than those in Fornax and even Sculptor. 

The SAAO reports that all three of our targets in Ca- 
rina are semi-regular variables, but the periods are un- 
determined (J. Menzies, P. Whitelock, and M. W. Feast, 
2011, private communication). As with Leo I, the metal- 
hcities in Carina are distributed more narrowly than in 
Sculptor or Fornax: <[Fe/H]> = -1.73 ± 0.35, making 
it likely that these targets are even more metal-poor than 
those in Leo L The metallicity measurement by Abia et 
al. (2008) reinforces this point; they found [M/H] = -1.7 
for MCA C3I8. 

Thus the metalhcities of the carbon stars observed in 
Fornax and Sculptor are probably higher than previously 
proposed. While Matsuura et al. (2007) adopted [Fe/H] 
^ — 1.0 for the five Fornax stars they examined, we find 
that [Fe/H] appears to be more Magellanic in nature: 
0.3 to -0.8. Similarly, while Sloan et al. (2009) sug- 
gested that MAG 29, which is also associated with signif- 
icant quantities of carbon-rich dust, would have formed 
with [Fe/H] ~ —1.4, we have revised this metallicity up 
to ~-1.0. 

Figure 6 includes the P-L relation of Whitelock et al. 
(2009): Mboi = -3.3 log P (days) + 3.979. This relation 
appears on the plot as a narrow dotted curve, but for typ- 
ical data, it is ~0.4-0.5 magnitudes wide at any period. 
Most of our data fall within this range. The evolutionary 
tracks show why a significant spread in periods would be 
expected for stars of a given mass (and therefore age). 
For stars of mass > 1.5 Mq, once they reach pulsation 
periods of ~350-400 days, their period will continue to 
increase, but their luminosity will remain largely fixed. 
Thus the width of the P-L relation depends on how long 
AGB stars will survive once they begin pulsating in the 
fundamental mode. This width can lead to considerable 
uncertainty in any distances determined using P-L rela- 
tions for LPVs, and for this reason we did not include 

1* And [M/H] = -1.9 for ALW Car 7. 



TABLE 9 
Fitting wavelengths 



Feature 


A 


Blue continuum 


Red continuum 




(/im) 


(fim) 


(Atm) 


C2H2 abs. 


7.5 


6.08-6.77 


8.25-8.55 


SiC dust cm. 


11. 


9. .50-10.10 


12.50-12.90 


C2H2 abs. 


13.7 


12.80-13.40 


14.10-14.70^ 



* 14.10-14.17 for all but the two spectra with LL data, 
distances based on the periods of LPVs in § 2.1 above. 



4.3. Comparison samples 

In the analysis below, we compare the carbon stars 
we have observed in the four targeted dwarf spheroidal 

galaxies with similar infrared spectroscopy of samples in 
the Milky Way, LMC, and SMC. The Galactic spectra 
are from the atlas of spectra from the Short- Wavelength 
Spectrometer (SWS) on the Infrared Space Observatory 
(ISO) (Sloan et al. 2003). The sample includes the 37 
sources classified by Kracmcr ct al. (2002) as carbon 
stars, including nine stars observed multiple times. Us- 
ing the multiple observations, we have arrived at typi- 
cal variations over the pulsation cycle of the star in the 
strengths of the extracted features; these are plotted in 
the relevant figures to give an idea of the systematic un- 
certainties in the data. See the description by Sloan et 
al. (2006) for more detail on the Galactic sample. A pa- 
per in preparation will present the spectral properties of 
this sample in more detail. 

The Magellanic samples of carbon stars come from the 
following programs: 200, 1094, 3277, 3426, 3505, and 
3591. The publications arising from these programs are 
referenced in §1. These samples include a total of 72 
carbon stars in the LMC and 34 in the SMC. 

Some caution is required when comparing spectral data 
amongst the samples, as a variety of selection criteria bias 
them in different ways. For example. Programs 1094 and 
3591 focused on dusty sources in the LMC, thus selecting 
against optically thin dust shells, while Program 3426 
sampled the brightest infrared sources in the LMC, which 
again selected against optically thin dust shells. The 
impact of these particular criteria is readily apparent in 
the following figures. For this reason, we will attempt 
to compare the data in the different samples by plotting 
the variable of interest against a dependent variable. The 
effect of metallicity will reveal itself through changes in 
the dependency from one galaxy to the next. 

5. SPECTRAL ANALYSIS 
5.1. The Manchester Method 

The mid-infrared spectra of carbon stars are rich in 

emission and absorption features. Figure 7 ilhistratcs 
how we sort these out for one spectrum, using the Manch- 
ester Method, which was first developed for IRS spectra 
of carbon stars in the Magellanic Clouds (Sloan et al. 
2006; Zijlstra et al. 2006). The [6.4] -[9.3] color samples 
the spectra at 6.25-6.55 and 9.10 9.50 /xm. These two 
wavelength intervals are relatively free of emission or ab- 
sorption features, and as detailed in § 5.2 below, can be 
used to estimate the total amount of dust emission in 
the spectrum. For ease of reference, we will refer to this 
quantity as "dust content" in the following discussion. To 
measure the relative strengths of the acetylene bands and 
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TABLE 10 
Spectral dust properties 



Source 


[6.4.] -[9.3] 


Mass- loss rates (Mq /yr) 


"SiC" feature strength 


name 


(mags.) 


total (logM)^ dust (logt))^ 


Ac (/irn) feature/continuum 



MAG 29 


0.432 


± 


0.008 


-5, 


.91 


-8, 


.21 


10, 


,87 


± 


0.08 


0.028 


± 


0.005 


Scl V78 V544 


-0.147 


± 


0.016 


















0.061 


± 


0.034 


Fnr RW 2 


0.333 




0.015 


—6 


.07 


—8 


.37 


11. 


,51 


-|- 


0.16 


0.084 




0.013 


For BTH 13-23 


0.716 


± 


0.006 


-5 


.45 


-7 


.75 


11, 


,29 


± 


0.07 


0.133 


± 


0.006 


For BTH 12-4 


0.602 


± 


0.005 


-5 


.64 


-7 


.94 


11, 


,27 


d= 


0.12 


0.045 


± 


0.004 


For BTH 3-129 


0.592 


± 


0.006 


-5 


.65 


-7 


.95 


11, 


,23 


± 


0.04 


0.173 


d= 


0.005 


For DK 18 


0.152 


± 


0.021 


-6 


.36 


-8 


.66 


11, 


,48 


± 


0.28 


0.120 


± 


0.024 


For DK 52 


0.198 




0.031 


-6 


.28 


-8 


.58 


11, 


,50 


± 


0.43 


0.093 


± 


0.031 


For DI 2 


-0.153 


± 


0.068 










11, 


.34 


± 


0.15 


0.317 


± 


0.053 


For WEL CIO 


-0.046 


± 


0.022 










11, 


.80 


± 


0.33 


0.051 


± 


0.015 


For BW 69 


0.232 


± 


0.023 


-6 


.23 


-8, 


.53 










0.011 


± 


0.026 


For BW 75 


0.156 


± 


0.019 


-6 


.35 


-8, 


.65 










0.032 


± 


0.018 


For BW 83 


0.163 


± 


0.032 


-6 


.34 


-8 


.64 


11, 


.93 


± 


0.26 


0.041 


± 


0.014 


ALW Car 2 


-0.025 


± 


0.028 


















0.012 


± 


0.047 


Car MCA C3 


-0.157 




0.018 


















-0.008 


± 


0.039 


Car MCA C5 


-0.031 




0.031 


















0.088 


± 


0.056 


Leo I MFT C 


0.500 


± 


0.015 


-5 


.80 


-8 


.10 


11, 


,24 


± 


0.45 


0.155 


± 


0.043 


Leo I MFT A 


0.069 


± 


0.016 


-6 


.49 


-8 


.79 










0.080 


± 


0.038 


Leo I MFT E 


0.226 




0.012 


-6, 


.24 


-8 


.54 


11, 


,32 


± 


0.22 


0.196 


± 


0.026 



Assuming a gas-to-dust ratio {ip) of 200. 
Assuming an outflow velocity {vout) of 10 km/s. 



TABLE 11 
Spectral acetylene properties 



Source 


7.5 /im C2H2 band 


13.7 /im C2H2 


band (Q branch) 


name 


Ac (/im) EW (jum) 


Ac (/im) 


EW (//m) 



MAG 29 


7, 


,49 


± 


0, 


,02 


0, 


,780 


± 


0, 


,012 


13.69 


± 


0.04 


0, 


.079 


± 


0, 


.007 


Scl V78 V544 


7, 


,31 


± 


0, 


,15 


0, 


,186 


± 


0, 


,013 








0, 


.144 


± 


0, 


,031 


For BW 2 


7, 


,34 


± 


0, 


,09 


0, 


,112 


± 


0, 


,010 








-0, 


,002 


± 


0, 


,029 


For BTH 13-23 


7, 


,43 


± 


0, 


,02 


0, 


,238 


± 


0, 


,004 


13.68 


± 


0.03 


0, 


,042 




0, 


,004 


For BTH 12-4 


7, 


,52 


± 


0, 


,08 


0, 


,050 


± 


0, 


,004 


13.70 


± 


0.02 


0, 


,058 


± 


0, 


,003 


For BTH 3-129 


7, 


,42 


± 


0, 


,02 


0, 


,191 


± 


0, 


,005 


13.66 


± 


0.02 


0, 


,030 




0, 


,002 


For DK 18 


7, 


.46 


± 


0, 


,11 


0, 


,160 


± 


0, 


,011 








0, 


,040 


± 


0, 


,032 


For DK 52 












0, 


,117 


± 


0, 


,024 


13.80 


it 


0.08 


0, 


,157 


± 


0, 


,027 


For DI 2 


7, 


.24 


± 


0, 


,13 


0, 


,183 


± 


0, 


,035 








0, 


.122 


± 


0, 


.051 


For WEL CIO 


7, 


.47 


± 


0, 


,05 


0, 


,178 


± 


0, 


,008 








0, 


.095 




0, 


.013 


For BW 69 


7, 


.44 


± 


0, 


,15 


0, 


,149 


± 


0, 


,021 








0, 


.034 




0, 


.021 


For BW 75 


7, 


.52 


± 


0, 


,11 


0, 


,136 


± 


0, 


,009 








0, 


.019 


± 


0, 


.021 


For BW 83 


7, 


.48 


± 


0, 


,16 


0, 


,146 


d= 


0, 


,020 


13.66 


± 


0.09 


0, 


,130 


± 


0, 


,023 


ALW Car 2 












-0, 


,074 


± 


0, 


,023 








-0, 


,341 


± 


0, 


,034 


Car MCA C3 












0, 


,040 




0, 


,014 


13.70 


± 


0.09 


0, 


,150 


± 


0, 


,020 


Car MCA C5 












-0, 


,018 


± 


0, 


,018 








0, 


,127 


± 


0, 


,028 


Leo I MFT C 


7, 


,52 


± 


0, 


,05 


0, 


,315 


± 


0, 


,010 








-0, 


,008 




0, 


,020 


Leo I MFT A 


7, 


,43 


± 


0, 


,08 


0, 


,251 


± 


0, 


,016 








0, 


,083 


± 


0, 


,036 


Leo I MFT E 


7, 


,53 


± 


0, 


,14 


0, 


,075 


± 


0, 


,013 








0, 


,039 




0, 


,019 



SiC emission features in the spectra, we use line segments 
to estimate the continuum above or below the feature. 
Table 9 gives the wavelengths used to fit continua and 
measure the various features. For the molecular bands, 
wc report an equivalent width (EW); for the SiC dust 
emission, we report its total integrated strength, divided 
by the continuum underneath, as estimated by the fitted 
line segment. For all features, we also report a central 
wavelength, defined as the wavelength which bisects the 
integrated flux of the feature. The uncertainty in the 
central wavelength indicates the range possible given the 
uncertainty in the extracted strength. 



Table 10 presents the [6.4] — [9.3] color and the relative 
strength of the SiC dust emission feature at ~11.3 ^m. 
The SiC feature is covered below in § 5.3. Table 11 
presents the equivalent widths and central wavelengths 
of the absorption bands from acetylene gas, as described 
in § 5.4. 

5.2. Dust content and pulsation period 

The [6.4] — [9.3] color measures the total emission from 
the star plus shell in two wavelength ranges between the 
various absorption and emission features. Amorphous 
carbon, which dominates the dust around carbon stars 
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Fig. 7. — An example of the carbon-rich dust and molecular fea- 
tures in the spectrum of Leo I MFT C. We use the Manchester 
Method to measure the [6.4] — [9.3] color, the strength of the SiC 
dust emission feature at 11.3 fira (or other dust features in the 
vicinity), and the acetylene absorption bands at 7.5 and 13.7 fim. 
Note that the latter is just the Q branch of a much broader fea- 
ture. The method estimates the strength of the SiC and C2H2 fea- 
tures line segments, with the continua estimated over the ranges 
where the lines are thin and the features measured where the lines 
are thick. Table 10 gives the [6.4] — [9.3] colors and SiC feature 
strengths, while Table 11 gives the equivalent widths for the acety- 
lene bands. 
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Fig. 8. — The relation between [6.4] — [9.3] color and pulsation 
period for our sample of carbon stars in Local Group dwarfs, com- 
pared to samples from the Galaxy, LMC, and SMC. The uncer- 
tainties in [6.4] — [9.3] color are generally smaller than the plotting 
symbols, and it is reasonable to assume that the uncertainty in 
period is about 10%. The dashed curve is a line fitted to all of the 
stars with periods in the Milky Way, LMC, SMC, and Fornax. All 
four galaxies show a similar dependence of total dust content on 
pulsation period, despite their differences in metallicity. However, 
all five stars in Sculptor and Leo I lie below the line at a statisti- 
cally significant level, indicating that the impact of metallicity is 
revealing itself for the most metal-poor stars sampled. 



(see JVIartin & Rogers 1987), has no resonances in the 
infrared, so that the apparent "continuum" in a spec- 
trum is actuahy a combination of star plus dust. The 
[6.4] — [9.3] color provides a means of measuring the rel- 
ative combinations of these two components and thus 
serves as a proxy for total dust content. 

Groenewegen et al. (2007, their Fig. 7) calibrated the 
[6.4] — [9.3] color to total mass-loss rate by applying ra- 
diative transfer models to a large sample of IRS spectra 
of JVIagellanic carbon stars. They found a linear relation- 
ship between the log of the total mass-loss rate and the 



[6.4] — [9.3] color. Because their models used the same 
gas-to-dust ratio {ip = 200) and the same outflow veloc- 
ity {vout = 10 km/s) for all stars, their calibration of the 
[6.4] — [9.3] color actually ties it directly to the total mass 
of warm dust contributing in the 6-10 fim spectral re- 
gion. Dividing M by the gas-to-dust ratio gives the dust 
production rate (£), DPR, or dust JVILR), and dividing 
D by the outflow velocity gives a quantity proportional 
to the optical depth of the radiating dust, which we are 
calling the dust content. 
Sloan et al. (2008) presented the dust-color relation in 

terms of log D vs. [6.4] — [9.3]. Adding a term to account 
for the outflow velocity and correcting their typographi- 
cal error. 



logZ?(^)=lo, 



Vout 



lOkm/s 
X { - 8.9+ 1.6 ([6.4] - [9.3]) } 



(1) 



This equation makes no assumptions about outflow ve- 
locity, and it is free of any dependence on gas-to-dust 
ratio. We will return to the question of how these quan- 
tities vary with metallicity in § 6.5 below. 

Table 10 gives the [6.4] — [9.3] color for each source. In 
order to translate these colors into more familiar quanti- 
ties, Table 10 also provides D using the above equation 
and an assumed outflow velocity of 10 km/s and M as- 
suming ijj = 200. It is important to remember, though, 
that the [6.4] — [9.3] color actually measures the dust con- 
tent. 

Sloan et al. (2008) compared the dust content in carbon 
stars in Magellanic and Galactic samples by plotting the 
[6.4] — [9.3] color as a function of the pulsation period of 
the star (their Fig. 29). They found that the amount 
of dust as measured by the [6.4] — [9.3] color increases 
with pulsation period. The scatter in [6.4] — [9.3] color 
at a given period is substantial. Within this envelope, 
no dependency on metallicity was apparent between the 
samples from the Galaxy, LMC or SMC. 

Earlier publications of data from our Local Group sam- 
ple did not have the benefit of the periods determined 
from the SAAO, but we are now in a position to com- 
pare our Local Group carbon stars directly to the other 
samples. Figure 8 includes data for carbon stars in For- 
nax, Sculptor, and Leo I, along with the comparison 
sample from the Galaxy, LMC, and SMC. The over- 
all dependency of dust content with pulsation period 
is unchanged. The new Fornax data appear to follow 
the same dependency, although the period coverage is 
smaller. This similarity is consistent with our revised 
metallicity. The figure includes a line fitted to all of the 
Galactic, Magellanic, and Fornax data: 



[6.4] - [9.3] = -0.227 + 0.00169 P (days). 



(2) 



Interestingly, the five spectra from Sculptor and Leo I 
all lie below the fitted line in Figure 8. The mean differ- 
ence for these five is —0.185 magnitudes, with an uncer- 
tainty in the mean of 0.047 magnitudes. The compari- 
son data, including Fornax, show a standard deviation of 
0.209 magnitudes about the fitted line, and with a sam- 
ple of 121 objects with periods, the uncertainty in the 
mean is 0.019 magnitudes. Adding the uncertainties in 
quadrature gives an uncertainty in the difference between 
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the fitted line and tlie data from Sculptor and Leo I of 
0.051 magnitudes, making the difference between them 
and the other samples 3.6 a. 

In contrast, the mean difference between the Fornax 
data and the fitted line in Figure 8 is only 0.018 magni- 
tudes. Comparing this difference to a standard deviation 
of 0.241 and an uncertainty in the mean of 0.098 mag- 
nitudes shows that Fornax follows the Magellanic and 
Galactic samples. 

The difference for Sculptor and Leo I is statistically 
significant. The carbon stars in these two galaxies al- 
most certainly do not belong to the same population as 
Fornax, the SMC, the LMC, and the Galaxy {p value ~ 
0.00014). Nonetheless, we have only five deviant spectra, 
further verification with larger samples would be helpful. 
While we have revised previous estimates of the metal- 
licity of the carbon stars in Fornax upward to Magellanic 
values and the metallicity in Sculptor from —1.4 to —1.0, 
the Leo I sample restores the range of [Fe/H] sampled in 
this analysis down to ^ —1.35. We conclude that for 
the most metal-poor stars in our sample, the impact of 
the initial metallicity of the star on its future dust pro- 
duction as a carbon-rich AGB star is large enough to be 
noticeable in the infrared. 

5.3. Silicon carbide dust emission 
5.3.1. In the dwarf spheroidals 

The SiC feature at ~11.3 fim is generally weak in this 
sample. Only three sources show an unambiguous SiC 
feature, although two more have spectral structure con- 
sistent with it. In Table 10, seven features have a S/N ra- 
tio < 2.5. For these, the central wavelength is omitted be- 
cause it is meaningless; their continuum-subtracted spec- 
tra are essentially noise in this region. Figure 9 illustrates 
the extracted features in the vicinity of 11.3 11.5 /J,m 
for the remaining 12 sources which have an extracted 
strength with a S/N ratio > 2.5. 

The five probable or possible SiC features all have cen- 
tral wavelengths of 11.2-11.3 fim. Three of these are 
in Fornax (BTH 3-129, 12-4, and 13-23), and they were 
all discussed previously by Matsuura et al. (2007). The 
other two arc in Leo I. In both cases, the emission profile 
is not a perfect match for SiC, but given the noise, SiC 
is the most likely explanation. 

Fornax DI 2 has an emission feature which peaks at 
11.3 /xm and has a width and shape consistent with the 
out-of-planc C H solo bending mode in polycyclic aro- 
matic hydrocarbons (PAHs). The problem with this as- 
signment is that no other PAH features can be identified 
with any confidence. There is a hint of the 8.6 /Ltm fea- 
ture, but noise masks the 6.2 /zm feature, and acetylene 
absorption would hide the 7.5 7.9 /xm emission complex. 
No 12.7 fiia feature is apparent. 

Other features in Figure 9 can be placed in three 
groups, labelled "X", "Y", and "Z". 

MAG 29 has the sole "X" feature, an apparent emission 
feature centered at 10.9 /xm, but its strength amounts 
to only 3% of the continuum integrated over the same 
wavelength range. Given the strong molecular absorp- 
tion bands in this spectrum (§ 5.4), we suspect that the 
apparent peak at ~11 /im is simply the continuum be- 
tween molecular absorption bands. C3 at 10 /xm might 
explain the drop to the blue side of this "feature" (Zijl- 



stra et al. 2006). The drop to the red could be due to 
the broad wings of the C2II2 band centered at 13.7 /xm, 
which grow wider for higher gas temperatures. 

Three spectra show what we are calling the "Y" fea- 
ture, which peaks at ^11.5 /xm and is sharper and more 
symmetric than the SiC feature. These features are weak, 
and the spectra have poor S/N ratios, preventing any 
conclusive statements about their carrier. Nonetheless, 
it should be noted that graphite produces an emission 
feature in this spectral range (Draine & Lee 1984; Laor 
& Draine 1993). The absence of this feature in the spec- 
trum of IRC +10216 (Martin & Rogers 1987) and other 
Galactic carbon stars led to the currently favored model 
whore amorphous carbon, and not graphite, dominates 
the dust around carbon stars. The 11.5 /tm graphite fea- 
ture arises from C-C displacements between graphene 
sheets, and in laboratory data it is exceptionally nar- 
row due to the regular spacing and large extent of these 
sheets. Smaller sheets and irrcgiilarities in the lattice 
structure would fatten the feature, but whether or not 
it would have a shape like the observed "Y" feature is 
unknown. The presence of graphite in carbon-rich dust 
shells is an interesting possibility, but given the limited 
quality of the extracted features, any further speculation 
is unwarranted. 

Two spectra show possible emission features in this 
range which peak further to the red, in the 11.8-11.9 /tm 
region. We are calling this the "Z" feature. As with the 
"Y" feature, though, the data are noisy. One could argue 
that the feature in the spectrum of For WEL CIO is a 
very noisy example of SiC emission, while the feature in 
For BW 63 could be a noisy "Y" feature. 

5.3.2. Comparing the samples 

Figure 10 plots the strength of the SiC dust emission, 
normalized to the underlying continuum as a function 
of the [6.4] — [9.3] color. The SiC dust strength is inte- 
grated between 10.1 and 12.5 /tm and divided by the 
total "continuum" emission in the same interval, where 
"continuum" is the combination of emission from amor- 
phous carbon dust and the central star. 

Sloan et al. (2006) found tliat the relative strength of 
the SiC emission feature decreased as the metallicity of 
the sample decreased. In Figure 10, two different se- 
quences of relative SiC strength vs. total dust content 
are apparent, with a clear bimodality at colors of ~0.2- 
0.6. In this range. Galactic stars dominate the upper se- 
quence, with not one falling on the lower sequence. Once 
[6.4] — [9.3] exceeds ~0.6, the sequences begin to merge, 
although the difference between the Galactic and Magel- 
lanic samples is still evident. 

Generally, the carbon stars in the dwarf spheroidals 
follow the lower sequence in Figure 10. Fornax DI 2 
(with [6.4]-[9.3] = -0.15 and SiC/cont. = 0.32) is the 
most significant exception, showing an emission feature 
at 11.3 /im with strong contrast to the continuum but 
virtually no other dust. This emission feature looks 
more like PAHs than SiC (see Figure 9), but as discussed 
above, that identification is problematic and uncertain. 
The blue [6.4] — [9.3] color indicates that For DI 2 is vir- 
tually naked, making SiC dust unlikely. Another major 
exception is Leo I MFT E ([6.4] -[9.3] = 0.23, SiC/cont. 
= 0.20). In this case, the feature is most likely SiC and 
the exception appears to be real. 
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Fig. 9. — The extracted spectral features in the vicinity of the SiC dust emission feature for the 12 spectra where the extracted flux 
has a S/N ratio > 2.5, along with identifications of the features where possible. "X" refers to the apparent feature in the spectrum of 
MAG 29 centered at 10.9 fim and most likely due to the continuum between molecular absorption bands. "Y" labels features with a central 
wavelength at ~11.5 fira, and "Z" labels features centered to the red, ~11.8-11.9 ^m. 

While the apparent SiC emission in Leo I MFT C is 
somewhat noisy, its profile is consistent with SiC. The 
identifications of the features in the two Fornax spectra 
are firm, due to the strength and profiles of these fea- 
tures. We see no obvious characteristics to distinguish 
these sources from the others with weaker SiC features, 
much as Sloan et al. (2006) found when investigating the 
five SMC sources in the same region. None of these eight 
stand out in terms of pulsation period, luminosity, or 
any other identified property. Why they have stronger 
SiC features than other sources from the same galaxies 
remains unknown. 

The sources in Fornax show SiC strengths similar to 
the Magellanic sources, fully consistent with the metal- 
licities we believe they formed with. The general lack of 
dust in the three Carina sources places them in the lower 
left corner of Figure 10, where the two tracks converge, 
giving us little new insight. The Leo I sources, however, 
are a bit of a surprise, with one of the three on the upper 
sequence. 
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Fig. 10. — The strength of the SiC dust emission at ~11.3 fim, 
relative to the underlying continuum from star and amorphous 
carbon dust, plotted as a function of [6.4] — [9.3] color. The error 
bars labelled "Typical variation" indicate the expected changes 
in spectral properties over a pulsation cycle of the central star. 
The uncertainties are plotted for the program stars, but they are 
generally smaller than the plotting symbols. The Galactic carbon 
stars follow a different sequence than most of the carbon stars in 
the other galaxies. The one data point from Fornax with [6.4] — [9.3] 
= —0.15 and SiC/cont. = 0.32 is For DI 2, which appears to show 
emission from PAHs and not SiC at ~11.3 fim. 

Three carbon stars in dwarf spheroidals appear in 
Figure 10 with SiC/cont. > 0.13 and [6.4]-[9.3] > 
0.45. These objects are (left to right in Figure 10): 
Leo I MFT C, For BTH 3-129, and For BTH 13-23. 



5.4. Acetylene gas absorption 

Table 11 presents the equivalent widths of the acety- 
lene (C2H2) bands centered at 7.5 fim and 13.7 ^m. As 
with Table 10, we do not quote the central wavelength if 
the S/N ratio of the equivalent width is less than 2.5. In 
addition, we omit it if the uncertainty in the wavelength 
exceeds 0.2 /xm at 7.5 /zm or 0.10 /im at 13.7 fim. In the 
cases with no central wavelength in Table 11, the equiv- 
alent width of the feature plus the uncertainty could be 
considered as an upper limit. 
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Fig. 11. — The equivalent width of the absorption band from 
C2H2 at 7.5 fim as a function of the [6.4] — [9.3] color. The data 
point in the top center is MAG 29 in Sculptor. This spectrum has 
a 7.5 fj,m band more than twice as deep as any other measured. 
The general trend of stronger absorption bands with decreasing 
metallicity is most apparent for 0.5 < [6.4] — [9.3] < 1.0. 

The 7.5 fiia band arises from the P and R branches 
of the vl + transitions and often presents a double- 
troughed structure with a central peak at 7.5 /im. The 
sharp absorption feature at 13.7 fj,m arises from the Q 
branch of the transition, primarily the fundamental 
mode, but with some contribution from higher overtones. 
The P and R branches extend this feature to cover the 
~13-15 ^m range, but our lack of LL coverage prevents 
us from measuring the full band. Table 11 presents the 
equivalent width of just the Q branch at 13.7 /im. Mat- 
suura et al. (2006) found that HCN, which often produces 
bands in the immediate vicinity of the C2H2 bands in 
Galactic carbon stars, was absent in the spectra of car- 
bon stars in the LJVIC. This conclusion also applies to 
carbon stars in the SIVIC (Sloan et al. 2006). No evi- 
dence of HCN appears in the Local Group spectra pre- 
sented here, but we note that limited S/N ratios and 
wavelength coverage at 14 fim prevent us from drawing 
any firm conclusions. 

The wavelengths in Table 9 were used for all features 
except one, the 7.5 /um band in JVIAG 29. That band 
is so strong that the continuum wavelengths fall in the 
wings of the absorption, forcing us to shift them to 6.14- 
6.44 /im on the blue side and 9.12-9.52 /im to the red. 
It is also possible that the extension of the absorption to 
longer wavelengths is due to another molecule, but that 
molecule would probably be another hydrocarbon and 
thus trace the same molecular hydrocarbon mix as the 
acetylene feature. 

We have detected a 7.5 fim band in most of our spectra. 
The exceptions include all three carbon stars in Carina 
as well as For DK 52. All of the spectra with detections 
show a clear minimum in the wavelength range covered 
by the band, although the noise can be considerable in 
some cases, due most likely to mismatches in the indi- 
vidual nod spectra and non-gaussian contributions which 
become a problem at these low flux levels. 

The narrow 13.7 /xm band is more difficult to detect. In 
addition to the lower signals at these wavelengths, only 
two of the spectra have LL data, leaving only three data 
points in SL for fitting the red continuum for the rest 



TABLE 12 

Mean SiC and acetylene strengths 



Galaxy 


Number 


<SiC/cont.>'' 


<7.5 Atm EW>=' 


Milky Way 


8 


0.29 ± 0.07 


0.06 ± 0.04 


LMC 


23 


0.11 ± 0.07 


0.08 ± 0.05 


SMC 


12 


0.08 ± 0.06 


0.12 ± 0.05 


For dSph 


3 


0.12 ± 0.06 


0.16 ± 0.10 



In the range 0.5 < [6.4] -[9.3] < 0.8. 



of the spectra, since we do not use SL past 14.17 /im. 
Nonetheless, we still detect a 13.7 /im band in seven of 
the 19 spectra, an impressive result given the sub-mJy 
strengths of most of the spectra at 14 fj,m. 

All of our sources in Sculptor, Fornax, and Carina have 
been confirmed as carbon stars in the literature. Held 
et al. (2010) describe two of the three targets in Leo I, 
JVIFT C and A, as probable carbon stars. IVIFT C shows 
a convincingly strong 7.5 /im absorption band, and the 
11.2 /im emission feature, while noisy, shows a shape con- 
sistent with SiC, allowing us to confirm its carbon-rich 
nature. IVEFT A also has a strong 7.5 /im band, although 
it is noisy. Given the lack of any obvious features at 11.3 
and 13.7 /im, this object remains an unconfirmed carbon 
star. The chemistry of IVIFT E is less clear in the litera- 
ture. The IRS spectrum shows what is best described as 
a weak, but noisy, 7.5 fim acetylene band and an emis- 
sion feature at 11.3 /im consistent with SiC, which leads 
us to treat it as a probable, but unconfirmed carbon star. 

We will focus on the 7.5 fj,m absorption band from 
C2H2 due to its better S/N ratio. Figure 11 plots its 
equivalent width vs. [6.4] — [9.3] color, allowing us to com- 
pare the molecular band strength at similar overall dust 
contents. The differences between the samples are most 
apparent in the color range 0.5-1.0, with a clear trend in 
increasing band strength from the Galaxy to the LIVIC to 
the SIVIC. Again, the three equivalent widths from For- 
nax are generally consistent with the IVIagellanic sample. 
Two of the three are among the strongest bands in this 
color range. 

The data from the other dwarf galaxies are outside this 
color range, but a couple of comments are in order. First, 
JVIAG 29 in Scl has an equivalent width more than twice 
as strong as anything else in the sample. IVIenzies et 
al. (2011) noted that the IRS spectrum was obtained at 
maximum luminosity, which might account for the strong 
acetylene absorption. Even if this band varied by a factor 
of two over a pulsation cycle, it would still be stronger 
when at its minimum than in any other spectrum con- 
sidered here. Second, two of the three stars in Leo I are 
among the strongest absorbers at 7.5 fj,m. The strong 
acetylene absorption in IVIAG 29 and in Leo I is consis- 
tent with the metal-poor nature of Sculptor and Leo I. 

5.5. Metallicity diagnostics 

Table 12 summarizes the mean and standard devia- 
tion of the strength of the SiC emission (normalized to 
the continuum) and the equivalent width of the 7.5 fim 
acetylene absorption band. Both of these quantities vary 
with [6.4] — [9.3] color, and the various samples consid- 
ered are incomplete at some colors. Consequently, we 
limit the comparison to the range 0.5 < [6.4] — [9.3] < 0.8. 
From the JVIilky Way to the LMC to the SIVIC, the trend 
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Fig. 12. — The equivalent width of the acetylene band at 7.5 fim 
plotted as a function of the SiC/continuum emission ratio. MAG 29 
is missing from this plot because it is so far above the rest of the 
data, with EW = 0.78 fim and SiC/cont. = 0.03. Lagadec et al. 
(2008) proposed this comparison as a metallicity diagnostic (their 
Fig. 11). 

of decreasing SiC strength and increasing C2H2 strength 
is clear. It is possible that the overlapping spreads in the 
data arise from a range of metallicities within the three 
galaxies. As described above (§ 2.2 and 4.2), we suspect 
that the metallicity of Fornax is closer to Magellanic than 
to the other dwarf spheroidal galaxies considered here. 
We have only three spectra in this color range, and they 
give ambiguous results. The SiC feature suggests a simi- 
larity to the LMC, but the spread in acetylene strengths 
limits our conclusions. 

Lagadec et al. (2008) proposed plotting the two quan- 
tities in Table 12 against each other to diagnose metallic- 
ity. Figure 12 follows that lead. While the sample from 
each galaxy shows considerable scatter, the gradient is 
unmistakable, with the Galaxy dominating in the lower 
right, the LMC in the middle, and the SMC in the upper 
left. The Fornax spectra are distributed much like the 
SMC spectra. As explained above, the one discrepant 
point in Fornax is For DI 2, which might have PAHs, 
not SiC, in its spectrum. Overall, Figure 12 reinforces 
our suspicions about Fornax and suggests that it and the 
SMC have similar metallicity distributions. 

The Sculptor data keep to the metal-poor side of the 
diagram, with MAG 29 literally off the charts (not plot- 
ted, with SiC/cont. = 0.03 and EW = 0.78 fim). This 
behavior is fully consistent with the narrow and metal- 
poor MDF for Sculptor. The three Carina spectra are 
clustered in the lower left, where metallicity is indeter- 
minant. The Leo I data are more of a puzzle. Two of the 
three spectra are in the upper left, where we would ex- 
pect them from the metallicity of their host galaxy, but 
Leo I MFT E lands squarely among the metal-rich spec- 
tra. This might indeed indicate a higher initial metal- 
licity, but it may also suggest that a bit of caution is 
warranted with these data. One alternative explanation 
is that this spectrum is affected by strong C3 absorption 
at 10 /.tm, which would push our continuum fit downward 
and enhance the apparent strength of the SiC emission. 

6. DISCUSSION 
6.1. Carbon dust content and metallicity 



The most surprising result of this study is that the car- 
bon stars in Sculptor and Leo I seem to have less dust 
around them than their counterparts in more metal-rich 
galaxies. Previous studies of carbon stars in the Galaxy, 
the LMC, and the SMC showed that the amount of cir- 
cumstellar dust, as measured by the [6.4] — [9.3] color, fol- 
lowed the same relationship with pulsation period of the 
star, despite their differences in metallicity. The six car- 
bon stars in Fornax with known pulsation periods also 
follow this relation. As explained in § 4.2, we believe 
that these stars have metallicities similar to those in the 
LMC and SMC. The two carbon stars in Sculptor, for 
which we estimate [Fe/H]'-^— 1.0, and the three carbon 
stars in Leo I, with [Fe/H]~— 1.35, are shifted downward 
on average about 0.19 magnitudes in Figure 8. This dif- 
ference corresponds to a factor of two in dust production 
(I), Equation 1). 

The scatter in Figure 8 is considerable. The standard 
deviation about the fitted line corresponds to a range in 
D of a factor of 2.2 (up or down). The cause of this 
wide envelope around the general proportionality of dust 
content to pulsation period is an open question. A more 
careful look at the Magellanic data rules out luminos- 
ity (Sloan et al., in preparation). Metallicity also does 
not work, even though we expect a fairly broad range of 
metallicities within each of the Milky Way, LMC, and 
SMC. If metallicity were responsible, it would lead to a 
measurable shift in the mean dust content between each 
galaxy, which we do not see. 

The data leave us to conclude that for [Fe/H] > —0.7, 
the amount of carbon-rich circumstellar dust shows no re- 
lation to the metallicity. The remainder of the discussion 
focuses on the implications. As stated in § 5.2 above, the 
carbon dust content is measured by the [6.4] — [9.3] color, 
with no assumptions about outflow velocity or gas-to- 
dust ratio. 

6.2. C/0 ratio and metallicity 

For a carbon star, the dust content depends primarily 
on the abundance of C, which it produces via the 3-q; 
sequence and dredges to its surface. The free carbon 
available for dust production depends on what remains 
after CO consumes the available oxygen, so that 



o o 



(3) 



In this equation and the ones that follow, "C" and "O" 
are the fractional abundances of carbon and oxygen, by 
number. Not all of the free carbon will form dust; some of 
it appears in our spectra as acetylene, and other carbon- 
bearing molecules are possible. Because our carbon stars 
are embedded in dust and their distances make them 
faint, they are difficult targets for the high- resolution op- 
tical spectroscopy usually used to determine C/O ratios. 
While waiting for that problem to be solved, we can still 
estimate what C/O ratios we should expect. We can 
write 

C Ci + SC fC\ f5C\ 



O 



O 



O 



O, 



where Ci is the initial carbon abundance, 5C is the 
change in abundance due to dredge-up, and we have as- 
sumed that the oxygen abundance does not change from 
its initial level. 
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Observations of the C/0 ratio in dwarf galaxies reveal 
a decrease in C/0 with decreasing rnetallicity (Garnett 
et al. 1995; Kobulnicky & Skillnian 1998, where 0/H 
serves as a proxy for rnetallicity). This is the initial C/0 
ratio prior to any dredge- up on the AGB. Data from stel- 
lar spectroscopy in the Galaxy compiled by Chiappini 
et al. (1999, 2003) show no significant trend in [d/Fe] 
with [Fe/H]. The scatter is substantial, ^0.2 dex, and the 
point is not without controversy, but we will assume no 
systematic relation between [Cj/Fe] and [Fe/H]. Similar 
spectroscopy by Melendez & Burbuy (2002, and refer- 
ences therein) shows that a steadily increasing [0/Fe] 
ratio in more metal-poor samples is the culprit behind 
the changing (C/0)i ratio. We estimate that in the range 
-1.5 < [Fe/H] < 0.0, [0/Fe] increases 0.25 dex for every 
decrease of 1 dex in [Fe/H]. The scatter about this trend 
is ~0.1 dex on average. Applying this trend to the above 
equation gives: 

(lOO-25 [Fc/H] ^ ^Q-0.75 [Fe/H] j _ (5^ 

The first term in the right-most parentheses is just the 
initial C/0 ratio, which decreases by a factor of 1.78 for 
every dex in [Fe/H]. 

Let us assume that the amormt of dredged up carbon 
(SC) does not vary significantly with rnetallicity. The 
arguments which follow will pursue the consequences and 
show that in fact this assumption is probably invalid. 

For the sake of argument, we will also assume that a 
typical C/0 ratio for a carbon star at solar metallicity 
is ~1.1 (following Lagadec & Zijlstra 2008, but see also 
Ohnaka et al. 2000). At solar mctalficity, the initial C/0 
ratio of a star would be 0.54 (Asplund ct al. 2005). To 
raise the C/0 ratio to 1.1, SC would have to be 0.56 Oq. 
Keeping this quantity fixed with metallicity leads to a 
C/0 ratio at analogous stages of AGB evolution of 1.4 
at [Fe/H] = -0.3 (i.e. in the LMC), 2.2 at [Fe/H] = -0.7 
(SMC), and 3.5 at [Fe/H] = -1.0 (Sculptor). 

Dividing Equation (2) by (C/O)© leads to a relation 
for the free carbon as a fraction of the total initial carbon 
at solar metallicity: 

("g-l) 1.85x100-^5 [Fe/H] _ (g) 

For a star of solar metallicity, Cfree = 0.19 Cq (in this 
particular case, (7© is the amount of carbon the star 
formed with). The corresponding values for the LMC, 
SMC, and Sculptor are 0.44, 0.68, and 0.81. It is in- 
structive to consider the example of the SMC. Despite 
the fact that stars in the SMC form with only 20% of 
the carbon in their Galactic counterparts, if they dredge 
up the same amount of carbon, they will have 3.6 times 
more carbon free to form dust. This is a direct result of 
the reduced oxygen abundance, which limits how miich 
carbon is sequestered as CO. All things being equal, we. 
would expect to see more carbon-rich dust at lower metal- 
licities, and yet, our measurements of the carbon dust 
content remain relatively constant with metallicity. We 
even see a possible dip at the metal-poor end of our sam- 
ple. 

6.3. Acetylene, the carbon budget, and mass-loss 
triggers 



Before giving up on our assumption that SC does not 
depend on metallicity, we should consider what might 
happen to this excess of free carbon at low metallicity. If 
the dust content stays flat or even falls, some other reser- 
voir must be absorbing the carbon. SiC represents only 
trace amounts, but the strengthening acetylene bands at 
lower metallicity could be a solution. 

Table 12 shows that the mean strength of the 7.5 /im 
acetylene band is twice as strong in the SMC than in 
the Milky Way, although the spread of the data about 
these means is considerable. If the 7.5 /xm acetylene band 
accurately traces the total carbon mass in gaseous hydro- 
carbon molecules, then it follows from Table 12 that hy- 
drocarbon gases sequester two times more carbon in the 
typical carbon star in the SMC compared to the Milky 
Way. 

While this is a significant amount of carbon in gaseous 
form, it is not enough to account for the factor of 3.6 more 

free carbon expected in the SMC. We know of no other 
likely reservoir for the carbon, leading us to conclude 
that the amount of dredged up carbon {6C) must in fact 
decrease at lower metallicity. Neither the rate at which 
the triple-a sequence produces carbon nor the dredge- 
up efficiency should decrease at lower metallicity. To 
reduce SC, it is necessary to terminate the dredge-up 
process progressively earlier on the AGB in more metal- 
poor stars. 

Models by Woitke (2006) show that the opacity of 
carbon-rich dust is sufficient to drive the mass-loss pro- 
cess. If a dredge-up event during a thermal pulse drives 
the C/0 ratio high enough, then the jump in free car- 
bon will lead to a pulse of dust formation which could 
quickly strip the envelope and bring the AGB evolution 
to a rapid end. Lagadec & Zijlstra (2008) proposed just 
such a scenario, although their definition of the free car- 
bon differs from ours. The key point is that the trigger is 
not the C/0 ratio, but the quantity of free carbon, mod- 
ified, if seeds are important to the condensation process, 
by the abundances of elements such as Ti and perhaps Si. 
A spectroscopic census of the Sgr dSph galaxy (McDon- 
ald et al., in preparation) finds that the typical lifetime 
of a star as a carbon star may be similar to the time 
between thermal pulses. Once a star develops a carbon- 
rich envelope, the formation of dust and the loss of mass 
accelerates quickly. 

6.4. Si, seeding, and abundance 

Van Loon et al. (2008) explained the increasing acety- 
lene band strengths in more metal-poor stars as a conse- 
quence of the lower abundances of elements like Ti. The 
lower abundances would result in fewer available seeds 
like TiC for the condensation of amorphous carbon dust, 
making the process less efficient. 

SiC could play a similar role as a condensation seed. 
Sloan ct al. (2006) suggested that the decreasing SiC 
emission at lower metallicity would arise naturally from 
decreasing Si abundances. Lagadec ct al. (2007) and 
Leisenring et al. (2008) developed this idea, explain- 
ing the different tracks of SiC emission strength with 
[6.4] — [9.3] color in Figure 10 in the different samples 
as an increasingly delayed condensation of SiC with de- 
creasing metallicity. Other seeds requiring heavy ele- 
ments should behave similarly with decreasing metallic- 
ity. Without these seeds at lower metallicities, the less 
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efficient formation of amorphous carbon is likely to result 
in a greater fraction of the free carbon tied up in simpler 
molecules like acetylene. 

Observations of Magellanic planetary nebulae (PNe) 
by Bernard-Salas et al. (2009) create one difficulty for 
the idea that the reduced SiC strength in metal-poor 
carbon stars arises from reduced Si abundances. They 
found that strong emission features centered at ^11 /im 
in several spectra, and they attributed these features to 
SiC. We suspect that the 11 /xm feature in their PN spec- 
tra does not arise from SiC, but instead from a different 
carrier. The PN features are triangular in shape and 
peak at ^11.1 /im, compared to the more rounded SiC 
features which peak at ^11.3 /im. Additionally, the PN 
features are usually associated with features at ~16 /zm. 
Both the 11 and 16 fim features are often seen in spectra 
showing the still unidentified 21 /xm feature (Kraemer et 
al. 2002, Sloan et al. in preparation). 

6.5. Gas-to-dust ratio and outflow velocity 

For evolved oxygen-rich stars, the gas-to-dust ratio 
should depend strongly on metallicity, because the abun- 
dances of the elements which form oxygen-rich dust scale 
with metallicity (e.g., van Loon 2000). The outflow ve- 
locity depends on the radiation pressure on the dust and 
its coupling to the gas due to collisions, and it should 
also vary with metallicity. 

The situation for carbon stars is different, because they 
produce the carbon which drives the mass-loss them- 
selves. The spectra in this paper show that if the dust 
content decreases with lower metallicity, the dependence 
is weak and noticeable only for the most metal-poor stars 
observed. Following the carbon budget leads us to the 
same conclusion as Lagadec & Zijlstra (2008), that the 
free carbon serves as a trigger for the superwind phase. 

The lack of a strong dependence of carbon-rich dust 
content with metallicity makes it difficult to see how the 
gas-to-dust ratio or the outflow velocity should depend 
on metallicity. If radiation pressure accelerates the dust, 
and collisions between dust grains and gas molecules ac- 
celerate the gas, then the amount of dust would deter- 
mine both outflow velocity and total mass-loss rate, and 
as we have shown, the amount of dust docs not vary 
strongly with metallicity. If pulsations drive the mass- 
loss process, then the mass- loss rate will depend primar- 
ily on the pulsation period and amplitude, which are 
relatively insensitive to metallicity. Again, neither the 
outflow velocity nor the gas-to-dust ratio should depend 
strongly on metallicity. 

A rigorous test of these predictions awaits observations 
from the Atacama Large Millimeter Array (ALMA) . Less 
sensitive telescopes can be used to study carbon stars in 
the Galactic Halo. Lagadec et al. (2010) observed a sam- 
ple of six carbon stars in the Halo and found that they 
have smaller outflow velocities than carbon stars in the 
Galactic disk. Lagadec et al. (2012) obtained infrared 
spectra of four of these carbon stars with the IRS, and 
they found evidence suggesting that halo Carbon stars 
also have high gas-to-dust ratios. Thus, both outflow ve- 
locity and gas-to-dust ratio may depend on metallicity 
for carbon stars, after all, but some caution is in order. 
First, the metallicity of the observed carbon stars is un- 
certain; they may be ejected from the Galactic disk. Sec- 
ond, if the outflow velocity of the dust is higher than the 



gas, then the inferred gas-to-dust ratios may not differ 
from the disk sample as much. The Halo carbon stars 
raise some interesting possibilities, but these initial re- 
sults require follow-up. 

The review of CO line observations in Galactic car- 
bon stars by Schoier (2007) shows a range of possible 
outflow velocities, from a few km/s to over 20 km/s, 
with higher outflow velocities generally associated with 
higher total mass-loss rates. The relation between out- 
flow velocity and total mass-loss rate is consistent with 
increasing amounts of circumstellar dust, which leads to 
more efficient acceleration of the circumstellar envelope. 
Metallicity is not needed to explain the effect, even if the 
range of possible outffow velocities does inject a degree 
of caution when using Equation (1) or interpreting the 
mass-loss rates presented in Table 10. 

7. SUMMARY 

The sample of carbon stars beyond the Magellanic 
Clouds observed by the IRS on Spitzer consists of 19 
stars in the Sculptor, Carina, Fornax, and Leo I dwarf 
spheroidal galaxies. The invidual targets were chosen 
based on their red NIR colors or their pulsation proper- 
ties. The bolometric magnitudes of the two targets in 
Sculptor, their pulsation periods, and reference to recent 
modelling work lead \is to revise their estimated metal- 
licity ([Fe/H]) up to ~— 1.0. Similarly, we suggest that 
the carbon stars in Fornax have metallicities more similar 
to the LMC and SMC than to the other dwarfs studied. 
The mean metallicity in Leo I and Carina are ~ —1.4 
and —1.7, respectively. 

All carbon stars follow a general relation of increasing 
dust content as the pulsation period increases. Metallic- 
ity does not show any influence down to [Fe/H]~— 0.7, so 
that despite their differences in initial metallicity, carbon 
stars in the Galaxy, LMC, SMC, and Fornax with simi- 
lar pulsation periods produce similar quantities of dust. 
The five stars observed in Sculptor and Leo I produce 
less dust, at a significance level of 3.6 a. The carbon 
stars in Carina show no dust, but lack of pulsation peri- 
ods prevent us from comparing them to the stars in the 
other galaxies. 

The new carbon stars extend the previously detected 
trends of decreasing emission from SiC dust and increas- 
ing absorption from acetylene gas as the metallicity falls 
to even lower metallicities. The carbon stars in Fornax 
show SiC features and acetylene bands consistent with 
our revised metallicity. 

If the quantity of freshly produced carbon from the 
stellar interior to the photosphere does not depend on 
metallicity, the lower initial abundances of oxygen in 
more metal-poor stars would lead to substantially greater 
quantities of free carbon. Yet, as our observations have 
shown, they do not produce substantially larger amounts 
of amorphous carbon dust. While more metal-poor stars 
do produce more acetylene gas, it is not enough to ac- 
count for the expected increases of free carbon. We con- 
clude that some process must truncate the dredge-up ear- 
lier in more metal-poor samples, quite likely a superwind 
triggered by a free-carbon threshold which strips the star 
of its envelope and ends its evolution on the AGB. 

Because the amount of circumstellar carbon-rich dust 
observed does not depend strongly on the metallicity, we 
should not expect any strong dependencies of outflow 
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velocity or gas-to-dust ratio on metallicity as well. 
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